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FEEDBACK BOUNDARY STABILIZATION TO TRAJECTORIES FOR 3D 

NAVIER-STOKES EQUATIONS 

SERGIO S. RODRIGUES* 


Abstract. Given a nonstationary trajectory of the Navier-Stokes system, a finite-dimensional feed¬ 
back boundary controller stabilizing locally the system to the given trajectory is derived. Moreover 
the controller is supported in a given open subset of the boundary of the domain containing the fluid. 

In a first step a controller is derived that stabilizes the linear Oseen-Stokes system “around the 
given trajectory” to zero; for that a corollary of a suitable truncated boundary observability inequality, 
the regularizing property for the system, and some standard techniques of the optimal control theory 
are used. Then it is shown that the same controller also stabilizes, locally, the Navier-Stokes system 
to the given trajectory. 
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1. Introduction 


Let H C be a connected open bounded subset located locally on one side of its smooth boundary 
r = and let / C M be a nonempty open interval. The Navier-Stokes system, in / x 11, controlled 
through the boundary reads 


dtu + {u ■'V)u — uAu + Vpu + h = 0, divrt = 0, ^i|r = 7 + C (1) 

where C is a control taking values in a suitable subspace of square-integrable functions in T whose 

support in x is contained in a given open subset Tc C T. Furthermore, as usual, u = (tti, rt 2 , 1 x 3 ) 

and pu-, defined for (t, xi, X 2 , X 3 ) G I x Q, are the unknown velocity field and pressure of the 

fluid, u > 0 is the viscosity, the operators V and A are respectively the well known gradient and 
Laplacian in the space variables (xi, X2, X3), {u ■ V)x stands for (tt • Vxi, u ■ Vx 2 , u ■ VU 3 ), divtt := 
dxiUi + < 9 x 2 '*^2 + 9x^U3, and h and 7 are fixed external forces. 

Suppose we are given a targeted (reference, desired) solution u{t) = u{t, x) of 0 with I = (0, + 00 ) 
and C = 0. If ii is stationary, u{t) = u{0, x) = uq, then the problem of stabilization to uq is now 
quite well understood. Namely, it was proven that, for any initial function uq sufficiently close to uq 
one can find a square integrable control C £ 9i^((0, + 00 ), L^(r, M^)), such that the corresponding 
solution u{t), supplemented with the initial condition 

u{0,x) = uo{x) ( 2 ) 


is defined on [0, + 00 ) and u{t) goes to uq exponentially as time t goes to + 00 ; we refer the reader to 
the works |Bad09( IBarl2t IBL12t IBTllt IFurOlt IFurOdt |Ray06| |Ray07[ IRTIO] . 

Again for a stationary targeted trajectory u = uq, the analogous result hold also in the case of an 
internal control under Dirichlet boundary conditions: 


dtu + {u ■'V)u — uAu ^Vpu + h + rj = 0, divM = 0, tt|p = 0. (3) 

where now 7 is a control supported in a given open subset w O Q. For details we refer to |Bar03L 
IBLT061IBTU4] . 

Here, we are particularly interested in the case where the targeted trajectory u is nonstationary 
(i.e., u = u{t) depends on time), a situation that often can occur in real world applications, as in 
the case suitable (say non-gradient) external forces {h and 7 ) depend on time. Also, since they are 
important and often required in applications, we look for controls obeying some general constraints 
like to be given in feedback form, finite-dimensional, and supported in a given (small) open subset. 

In [BRSllj . an internal stabilizing finite-dimensional feedback controller was found for the case 
of nonstationary targeted solutions. Then, one question arises: can we find a similar boundary 
controller? The methods used in the particular case of a stationary targeted solution, use some 
(spectral-like) properties of the (time-independent) Oseen-Stokes operator u <—)■ vAu — B{uo)u + Vpu 
and/or of its “adjoint” v <—)• i'Av — B*{uq)v + 'Vp^, which seem to give us no hint for the nonstationary 
case. Here B{uq)v := {uq ■ V)v + {v ■ V)uo and B*{uo) is the formal adjoint of B{uo). 

Also the constraints on the boundary control, imply that the procedure in [FI99] . that allow to 
derive suitable boundary results from internal ones is (or may be) no longer sufficient to derive the 
wanted boundary stabilization result. 

Departing from an exact controllability result in |Rodl4) , suitable truncated boundary observabil¬ 
ity inequalities have been derived for the (linear) Oseen-Stokes system in |Rodl5| . These results 
will enable us to follow the procedure in [BRSllj in order to construct a boundary stabilizing finite¬ 
dimensional controller to a given nonstationary targeted solution. To prove that the control can 
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be taken in feedback form, and to find the feedback rule, we will need to overcome some technical 
regular!ty/compatibility issues. 

More precisely, we will find the feedback rule by considering, at a first step, the linear Oseen-Stokes 
system 


dtv + B{u)v — uAv + Vpy = 0 , (4a) 

divu = 0, u|p = Cj (4b) 


and we look for the control C that minimizes a suitable cost J{v, C); by the Karush-Kuhn-Tucker 
Theorem and the dynamical programming principle we will conclude that the control is given in 
feedback form. Though this is a standard procedure, also used in |BRSllj . in the boundary control 
case we meet some nontrivial regularity/compatibility issues. Roughly speaking in the internal case 
to have the feedback rule at time t we need to know q{t) G L^(0, M^) where is a suitable Lagrange 
multiplier, associated with the constraint ( [da] ), that solves a system “adjoint” to Q, while in the 
boundary case we will need to know n-Vg(t) on the boundary T, where n is the unit outward 

normal vector to the boundary T and Pq(t) is a function depending on q{t). Thus, we will need more 
regularity for q. One possible way to get more regularity for q is to consider a cost functional J(u, C) 
that penalizes u in a larger (i.e., less regular) space, but we need to keep enough regularity for the 
optimal solution v to give a meaning to v{t), because we recall we want the control C in feedback 
form, that is, we want C, in (4b), as a function of v{t): ((t) = K{t, v{t)). We shall guarantee enough 
regularity for both q and v by considering an appropriate cost functional and an appropriate auxiliary 
extended system (cf. [BadOQ] 1. 

We shall prove the following Theorem, whose exact formulation is given in Section 
Main Theorem: Let {u,pu) be a global smooth solution for problem ([^, with C = 0 and t G Mq = 
( 0 , +oo), such that 


|ii|L°°(]RoxO,]R3) + sup |5lR|L2((r,T+l),L'^(0,R3)) + SUp | V'u|2,2((,r, t+ 1), L3(n, R9)) < R 

re[0,+oo) Te[0,+oo) 

where ii > 0 and o' > | are constants. 

Then for any A > 0 and any open subset Tc C T = dLl there are an integer M = M{R, A, Tc) > 1, 
an M-dimensional space £m C {/ G (^^(r,]^^) | /lr\rc ~ ® family of continuous linear 

operators t>0, from a suitable subset of into £m, such that the following assertions 

hold. 

(a) The function 1 1 —)• is continuous in the weak operator topology, and its operator norm is 

bounded by a constant depending only on R, A, and Tc. 

(b) For any divergence free function uq G such that the difference uq — n(0) is suffi¬ 
ciently small in the H^{£l,M.^)-norm and (uq — 'u(0))|p G £m, problem with ({t) = 

{uo - h( 0 ))|r + /o kV {u — u{r))dr has a unique global solution {u,pu), which satisfies the 
inequality 

|tt(t) — JJ3) < de 1^0 ~'n(0)liyi(o^]R3), t > 0. 

At this point we should say that this Theorem remains true for the two-dimensional (2D) case. 
Though we will focus on the three-dimensional (3D) case, the procedure is still valid for the two- 
dimensional one. 
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The fact that the control appears in integral form is meaningful from the physical point of view; 
indeed, since u is the velocity of the fluid then, roughly, the integral form means that we are accel¬ 
erating (or forcing) the fluid particles through the boundary. From the physical and practical point 
of view this is more natural than instantaneously imposing the velocity of the boundary particles. 

Though the integral feedback form of the controller, we will also show that the control is defined 
pointwise in time, that is, the control C(t) at time t > 0 depends only on u{t) — u{t), and not on the 
trace (uq — ri(0))|p as the integral feedback form could suggest. 

The feedback control we are going to construct will have both tangent and normal components. 
In some particular cases, in the case the targeted trajectory u is stationary, the stabilization of the 
Navier-Stokes system by normal boundary controls is proven in [BLKOH IBarOT^ IMunl2bl IMunl2a1 
IVK05| . In the general case, stabilization to a stationary solution in the 2D case has been achieved 
in [Bar la, under some general conditions, by means of oblique controls. This oblique stabilization 
result also holds in the 3D case for the linear Oseen-Stokes system. We would like to refer also to the 
work [Barl3j where the idea in [BarI2j is used for boundary stabilization to a stationary solution of 
parabolic equations, and leads to a simple algorithm to construct the stabilizing controller. Finally for 
the stabilization to a stationary trajectory by means of tangential controls we refer to |BLT06LIBLI2j . 

The rest of the paper is organized as follows. In Section we introduce some functional spaces 
arising in the theory of the Navier-Stokes equations and recall some well-known facts. Sections 3 
and 1^ are devoted to studying the linearized problem, that is, the Oseen-Stokes system; in Section 3 
we prove the existence of a stabilizing control and in Sectionwe prove that the control can be taken 
in feedback form. Finally in Section we establish the main result of the paper on local exponential 
stabilization of the Navier-Stokes system. The Appendix gathers a few more remarks concerning 
some points in the main text. 

Notation. We write M, Z, and N for the sets of real, integer, and nonnegative integer numbers, 
respectively, and we define Ma := (a, -|-oo) for all a G M, and No := N \ {0}. We denote by D C 
a bounded domain with a smooth boundary F = dCl. Given a vector function v: {t,xi, X2, x^) i— 
v{t,xi, X2, xs) G M^, k G No, defined in an open subset of M x D, its partial time derivative ^ will 
be denoted by dtv. Also the spatial partial derivatives ^ will be denoted by dxiV. 

Given an open interval I C M, then we write W{I, X, Y) := {/ G X) \ dtf G T^(/, T)}, 

where the derivative dtf is taken in the sense of distributions. This space is endowed with the 
natural norm \f\w{i,x,Y) ■= (l/li2(/^x) + X = Y we write H^{I, X) := 

W{I, X, X). Again, if X and Y are endowed with a scalar product, then also W{I, X, Y) is. The 
space of continuous linear mappings from X into Y will be denoted by C{X Y). In the case 
X = Y we write simply C{X). 

C*[ai,...,afc] denotes a function of nonnegative variables aj that increases in each of its arguments. 

C, Ci, z = I, 2,..., stand for unessential positive constants. 


2. Preliminaries 

2.1. Functional spaces. Let D C be a connected bounded domain of class C°° located locally 
on one side of its boundary F = (9D, with Jp dF < -|-oo. 
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We recall some spaces appearing in the study of the system 0 (cf. |Rodl4[|Rodl5| l. We start by 
the Lebesgue and Sobolev subspaces 

M^) := {u G | divu = 0 in fl}, 1 < r < + 00 , 

M^) := {u G | divu = 0 in Q}, s > 0 . 

The incompressibility condition allows us to define the trace of n • n on the boundary T = di}, 
where n is the unit outward normal vector to the boundary T, and then to write 

H := {u ^ I • n = 0 on T}, He ■= {u G | u • n = 0 on T \ Tc}, 

where Tc is an open subset of T. Some spaces of more regular vector fields we find throughout the 
paper are 

V :={u£ M^) I u = 0 on T}, W := {u G M^) | u = 0 on T \ 1^}, 

D{L) :=VnH‘^in,R^). (5) 

The spaces M^) are endowed with the scalar product inherited from H^{Q,,R^); the 

spaces H and He with that inherited from the spaces V and Ve with that inherited 

from H^{Q,R^); and D(L) with that inherited from H'^{Q,,R^). Notice that denoting by IT the 
orthogonal projection in onto H, it is well known that D(L) coincides with the domain 

{u G V\Lu € H} oi the Stokes operator L := —zzIIA. That is the reason for the notation. 

Next, fix a constant u > |. For any pair of real numbers a, b, with a < b, we introduce the 


Banach spaces and of the measurable vector functions u = {ui,U 2 ,U 3 ), 

defined 

in (a, b) x D, satisfying 




/ 

u\l 


(6) 

W(a, i))|Bt •= ( 

^ ^(a,i))|wk + L2((a,6),L3(n,R9 ))J <00- 


and also the Morrey-like spaces 




yywk _ ^ 


SBp.,-g[0, -l-oo) 1^ 1 (r, T-l-1) 1 WG’■^+l)l"'k ^ Too|', 

(7) 

IW* := < 

“ 

sup.,-g[Q^ Irt 1 .r+i)|st < Too J, 


endowed with the norms u yywk 

— sup 1 U 1 'yy(T, r + 1) 1 wk ; Und Iwlyyst . — SUp | U | ■yyfr, t + 1) |st ■ 



r>0 T>0 


Remark 2.1. The lower bound | for a is motivated from the results in [FCGIPOdl IRodlhj . 

We recall that, in |FGH02j . the set of traces ri|p at the boundary F of the elements u in the space 
lT((a, b), M^), M^)) is completely characterized, for each s > 5 , with s ^ 

Denoting that trace space by 6 ), F), we have that n 1 —>■ n|p is continuous; 

and, there is an extension Eg'. b), F) — > W{{a, b), M^), M^)), which is con¬ 

tinuous: 

{Esw)\r = w and < C' 2 |u;|G|^((a_b),r)- 
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We will use only the cases s = 1 and s 


GL((a, b), r) 


2. From [FGH021 Section 2.2] we know that 
6 ),F)©G*,,,((a, 6), F)n, 


with 


Gi((a, 6), F) 
Gi,av((a, b), F) 
Gliia, b), F) 
<av((a, b), F) 


L\{a, b), Hkr, TV)) n Fl((a, 6), i/-i(F, TF)) 
L\{a, b), Hl{r, M)) n b), R-^F, M)) 

L2((a, b), Hl{r, TF))nHl{{a, b), H%F, TF)) 
L2((a, b), HliF, M)) n b), H~)(F, M)) 


where the subscript “av” stands for the zero averaged condition, that is, Hl^{F, M) := {/ G 
H^(F, M) I Jp/dF = 0}, and TF stands for the tangent bundle of the manifold F, that is, ele¬ 
ments of TF are (tangent) vector fields in F. 

For technical reasons we relax a little the trace spaces: we define the superspace G®((a, b), F) of 
Gl^{{a, b), F) by just omitting the average constraint: 


G^((a, b), F) := G|((a, 6), F) © G^((a, 6), F)n, 


( 8 ) 


with (cf. |Rodl5( Section 2.1]) Gn((a, b), F) := L‘^{{a, b), H^{F, M)) n 6), H ^(F, M)) and 

Gliia, b), F) := L^{{a, b), hUf, R)) n H^ia, b), R-^F, M)). 


2.2. The control space. Let us write L^(n,M^) = R (B R'^, where R-^ = {V^ j ^ G R^{Fl, M)} 
denotes the orthogonal complement of R in L^(n,M^), and denote by IT the orthogonal projection 

Fi: L‘^{n,R^) ^ R (9) 

in L^(n,M^) onto R. For each positive integer N, we now define the iV-dimensional space Rn C R 
as follows: let {e^ j i G No} be an orthonormal basis in R formed by eigenfunctions of the Stokes 
operator L, whose domain is defined by ([^, and let 0 < ai < 0:2 < • • • be the corresponding 
eigenvalues, Let = aiSi, then put 

Rn ■= spanjej j i < N} C D(L) C R (10) 

and denote by FItv the orthogonal projection Flpj: R ^ R^ \n R onto Rn- 

Let O C F be a connected open subset of the boundary F, localized on one side of its boundary. 
We suppose that O is a C°°-smooth manifold, either boundaryless or with C'^-smooth boundary 
do. Let {vTj j i G No} be an orthonormal basis in L‘^{0, M) formed by the eigenfunctions of the 
Laplace-de Rham (or Laplace-Beltrami) operator Ao on the smooth manifold O, under Dirichlet 
boundary conditions, ©(p) = 0 for all p G dO. Analogously let {r* j i G No} be an orthonormal basis 
in L^{0, TO) formed by the vector fields that are eigenfunctions of Aq on TO, also under Dirichlet 
boundary conditions in the case dO 0, Ti(p) = 0 G TpF for all p G dO. It is known that vrj and Tj 
(i G No) are smooth. Let 0 < /3i < /32 < ..., and 0 < 71 < 72 < ... be the eigenvalues associated 
with the systems {vTj j i G No} and {rj j i G No}, respectively. 

We may write L?'{0, M^) as an orthogonal sum L‘^{0, M^) = L^(0, M)n © L?'{0, TO). Notice 
that {©n j i G No} is an orthonormal basis for L?‘{0, M)n = {/n j / G L‘^{0, M)}, and the system 
{vrjn j i G No} U {r* j i G Nq} is an orthonormal basis in the space L‘^{0, M^). 

Define, for each M G No, the space 

L\^{ 0 , M^) := spanlvTiii, © ] i G No, i < M} 


( 11 ) 
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and, denote by the orthogonal projection P^: L‘^{0, —)• L\j{0, in L‘^{0, M^) onto 

LliO, _ 

As in |Rodl5l Section 2.2], we suppose that the control region Fc and O satisfy 


Fc = supp(x) for 

some X G G^(F, M); and Fc C Fc C O C F. 

(12) 

Let us define the control space 



Jl 

^3 

C{t) = 5«;(t), and k G H^{Ro, M^^)} 

(13) 

with 

H : M 2 M ^ 

:= spanlxE^i^^TTin, xE^Tj 1 z G No,i < M} 

(14) 

[ I ] 

Gi 

— 'Eiii X^o + ^M+iT^ ; 


where E ®: L‘^{0,R) —)• L^(r, M) stands for the extension by zero outside the subset O, and 
P^j^: L‘^{0, R^) —> {xnl^}-*- stands for the orthogonal projection in L‘^{0, onto {ynl^p}-*- = 
{/ £ I (/, xn)L2('(p r3) = 0}: 


v 


and 


P^^v :=v- 


XTl)i2((p,K3) I 
ioX^dO 


Remark 2.2. Notice that the function ( = Sz satisfies the zero-average compatibility condition: 
fr CW • ndr = YlZi fo ^y-L(nn) • ynl^ dO = 0. 


In particular, observe that the controls, in Em-, are supported in [0, + 00 ) x Fc and take their values 
Q{t) in the finite-dimensional space Em, for each t £ [ 0 , -|-oo). 

Let us be given a constant A > 0 and two (fixed) regular enough functions h and 7 ; in addition we 
suppose that u is also regular enough and solves, in Mq x the Navier-Stokes system ([^ with C = 0) 
and a suitable pressure function = Pu- Given u(0) close enough to u(0), with (u(0) — m(0)) |p G Em, 
our goal is to find a (time-dependent) feedback linear controller v 1 —)■ G such that the 

solution of the problem 0 . with ( = (u( 0 ) — u( 0 ))|p-|-E Jq K^’'’(u — u(r)) dr, is defined for alH > 0 
and converges exponentially to u, with rate that is, 

\u{t) — 'h(^)|j:fi.^(Q JJ 3 ) < Ge \u(0) — fl( 0 )|p^i^^(Q JJ 3 ) for t P 0, 

where C is independent of n(0) — u{0) and time t. It will be clarified later in Sectionwhat we mean 
by “regular enough”, “close enough” and “solution”. 

Notice that seeking a solution of in the form u = u + v, formally we obtain the following 

equivalent problem for v: 


(15) 


dtv + B{u)v + {v ■'V)v — uAv pS/pv = 0, divr = 0, 

u|r = C, r'(O) = uo-u(O), 

with Pv = Pu— Pu, and B{u)v stands for {u ■ V)v -|- {v ■ V)n. We can see that it suffices to study the 
problem of stabilization of system (15) to the zero solution. We shall start by deriving the (global) 
stabilization of the Oseen-Stokes system, in Mq x fl: 


dtv + B{u)v — vAv+ Vp-u = 0, divr = 0, 

r|r = C, n( 0 ) = uq, 

to the zero solution; from which we shall derive the local result for 


(16) 
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2.3. Weak and strong solutions, and admissible initial conditions. We briefly recall some 
notions and results from |Rodl4l IRodl5| concerning the weak and strong solutions for the Oseen- 
Stokes system, in a bounded cylinder (a, b) x 0,, with a, b real numbers, 0 < a < 6 . 


dtv + B{u)v — uAv + Vpv + (7 = 0 , 

^;|r = C = Kr], 

Recall the extensions s G {1, 2}, in Section [2d| 


div V = 
v{a) = 


0 , 

Vq. 


(17) 


Definition 2.3. Given u G vq G M^), g G b), M^)), and ( G 

Gl^{{a, b), r); we say that v, in the space W((a, b), M^), M^)), is a weak solution 

for system if y := V — Ei( G +^((a, b), V, V') is a weak solution for 


dty + B{u)y - uAy + Vpy + / = 0, div y = 0, 

y\r = 2 /(«) = yo 


(18) 


with f = g + dtEi( + B{u)E\Q — I'AEiC, and yo = vq — Ei({a) G H. Here weak solution for (18) is 
understood in the classical sense as in |Lio69t Chapter 1, Sections 6.1 and 6.4], [Tem95[ Sections 2.4 
and 3.2], |Temfll( Chapter 3, Section 3]. 


Definition 2.4. Given u G Vo G M^), g G L‘^{{a, b), M^)), and also C £ 

GUia, b), r); we say that v, in the space W{{a, b), M^), L^(H, M^)), is a strong solution 

for system (©, if y := V — E 2 C G R^((a, b), D(L), H) is a strong solution for system ( [l^ with 
f = g + dtE 2 C + B{u)E 2 C — VAE 2 C,, and yo = vo — E 2 C{a) G V. Again, strong solution for (18) is 
understood in the classical sense as in |Tem95[ Section 2.4]. 


In the case our control C is in the space £m a natural question is: what are the admissible initial 
vector fields vo, if we want to guarantee the existence of a weak solution? Notice that, from (13), ( 
takes the form ( = Hk with k G H^(M.o, M^^). It is also not hard to check that the mapping H, in 
the definition of the control space (13), maps H^{{a, b), M^^) into Gav((a, ^)) T) C ^)) T) 

continuously, that is, our control is of the form as in © with K = 3. 

The set of admissible weak initial conditions for system (©, with C G £m, is given by . 4=1 = 
H + 77 = 1 , ^=1 •= EiEH^{{a, b), M^^)(o) = { 7 (a) | 7 = EiEg and 7 G b), M^^)}. 

Similarly, the set of admissible strong initial conditions for system ©I with ( G £ m : is given by 
Ae 2 := 1 ^ + 77 = 2 , with 77=2 - E 2 EH\{a, b), R'^^){a). 

Moreover 77=i, M=i, 77=2 71=2 are Hilbert spaces, with associated range norms 



:= inf ■ 


1 u = EiEg{a), g G 77^((a, h), ^ 


^Uhi 

;= inf ■ 

|(w^, z)\hxHs^ I 

u = w + z and {w, z) G H x 77=i | , 

(19) 

uIa- 

1-^=2 

:= inf • 

X 

[I] 

to 

u = w + z and {w, z) G V x 77=2 1 • 



Theorem 2.5. If u e g G L‘^{{a, b), M^))^ and g G H^Ha, b), 

with Vo — EiEg(a) G H, then there is a weak solution v in IT((a, b), M^), 77~^(H, M^)) 

for system ©, with C = Eg. Moreover v is unique and depends continuously on the given data 
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(^^ 0 , g, v): 


|2 


- ^[I“lw(-,i-)|wk] ^ l^'li2((a,6),//-l(n,R3)) + \v\m{{a,b),R 2 M) 

Theorem 2.6. If u e g G L‘^{{a, 6), and g G H^{{a, b), 

with vq — E 2 Eg{a) G V, then there is a strong solution v in W{{a, 6), M^), L‘^{Q, M^)) for 

system with C = Eg. Moreover V is unique and depends continuously on the given data (vq, g, g): 

l „,|2 


m((a, b),Hl^^{n, R3), L^{n, R3)) 


- ^[l«lw(<^,6)Nt] (l^0lHi;^(n,R3) + l5lL2((a.fe),L2(n,R3)) + l^lHi((a,fe),R2M) 

Remark 2.7. The weak solution given in Theorem 2.5 does not depend on the extension Ei. Also, 
the set of admissible weak initial conditions is independent of Ei (cf. |Rodl4l Rems. 3.2 and 3.4]). 
Analogously, the strong solution given in Theorem 2.6, and the set of admissible strong initial 
conditions are independent of E 2 . 

2.4. Smoothing property. The following Lemma will play a key role. 


Lemma 2.8. Let us he given u G g G L‘^{{a, b), L^(Q, M^)), vq G -4=^, and also g G 

H^{{a, b), M^^), with vq — EfEg^a) G H; then for the weak solution v of system 0 with = Eg, 
we have (• — a)v G W{{a, b), M^), L^(n, R^)), and 

\i' ~ ^)'^\wi{a,b),Hl (Q,R 3 ),L 2 (q^ 


< Ci 


div 

[l^lwKMIst] l^liR(a,6),L2(n,R3)) + l^l^i((a, 6),R2M) 


Proof. Since v solves ([l7|), it turns out that also w = {• — a)v does, with different data: 


dtw + B{u)w — u/S.w + V(- — a)py + (• — a)g — v = 0, 


w 


divtc = 0, 

w(a) = 0. 


Then, from Theorem 
C 


2.6 


|2 

\L^{{a,b),M{n, 


- a)g - n'2 


< \v\ 


= B.{--a)g, 

we can derive that the norm \w\w(^(a,b),Hl in,RS),M{Q,R 3 )) is bounded by 
L3{{a,b),L3(n,R3)) + !(■ “ ®)^lHi((a,fe),R2M))’ Lemma 
H-3{n, 


Wi(a,b),Hf(n, 


and Theorem 


2.5 


2.8 


follows from 

□ 


3. The Oseen-Stokes system: existence of a stabilizing control 

Using a controllability result from |Rodl5| . we shall construct a control, in the space Sm, expo¬ 
nentially stabilizing the linear system (16) to the zero solution. 

Definition 3.1. We say that vq is a weak, respectively strong, admissible initial condition for 
system (16), with f = Ek G £m, if it is a weak, respectively strong, admissible initial condition for 
the same system in (0, 1) x 11 with G EH^({0, 1), 

Let be the set of admissible initial conditions for controls ^ in £m, that is, Ae^ = 

H + EiEH^{{0, 1), 
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We introduce the mappings —)> and '■ —>■ defined by Q^V '■= 

(zi, Z 2 , •.., Z 2 m) with Zi = yi \i 1 < i < M, and Zi = 0ifM + l<i< 2M; and Qf^ := 1 r2m — ■ 

That is, QY is the orthogonal projection onto the first M coordinates, and QY the orthogonal 
projection onto the last M coordinates. 

In this Section we prove the following: 


Theorem 3.2. Let us he given A > 0 and u G W®*. Then there exists M = C'|'|^| ^ > 1 with 

the following property: for each {vq, k^) G Ae^ x ; there exists a “control” vector function 

u,\ _ i^O'j £ such that the weak solution v of system (16) in Mq x with 


K 

f satisfies the inequality 




—At 


|2 I I 0|2 


t > 0 . 


Moreover, for 0 < A < A, the mapping (uq, kY ^ is linear and satisfies: 


|e2 


(-^o, K. 




(A-A) 




+ Ik 


0|2 


3.1. Auxiliary results. Since the trace k(0), at time t = 0, is well defined for any given k G 
i/1 (Mo, we can easily obtain the explicit form of the spaces Aei and Asj, of all admissible 

weak and strong initial conditions, for controls in the space £m defined in (13). 

Writing v G Aei as v = u + {EiEk){0), where u e H and k G 1), R^^), from the fact 

that EiEk G bk((0, 1), M^), M^)) C ^([O, 1], L‘^{Ll, M^), and the continuity of the 

mapping u ^ u ■ n from M^) into ii“ 2 (r, M), we obtain that {v ■ n)n = limi^o((-£'i^K) • 

n)(i)n = limt \ o {{ EiEK { t )) ■ n)n = limt \ o {{ EK { t )) ■ n)n = ((Hk( 0)) • n)n. 

It follows that 


Alsi = {n G L2-y(n, M^) I {u ■ n)n|p = EQYz, for some z G . (20) 

Analogously, we can conclude that the set of strong admissible conditions is given by 

Ae 2 = {'U G //jjy(n, M^) I u|p = Ez, for some z G . (21) 

Remark 3.3. Notice that for a given u G L?- Ill, M^), we may set (u ■ to., ib) i , i , := 

(m, V 7 ^p'i/l)£, 2 (f 2 ]r), where T^p : i/s (T, M) —)■ M) is a continuous linear right inverse of the trace 

mapping T i—?■ Tip; in particular we have (7^pV’)lr ~ mapping u u ■ n|p is linear and 

continuous: hr • nl i < \u\TT\7if\ / i n. See. for examole. ITemOll Chanter 1. 

' 'R-2(r,R) — I in I i i£(iy2(r,R)^-n'i(f^,®)) ‘^ 

Section 1.3] and |LM72[ Chapter 1, Section 8.2]. Actually, writing (u • n, r/i) _i i := 

M 2 (^r, M), 2 (r, M) 

(u, V 77 p' 0 )^ 2 (f^ K) + (divu, '17.pV^)^2(Q K) we can see that the trace [v ■ n)|p G E{L\ —>■ //” 2 (r, M)) 
is well defined, with := {z G LfiVt,, M^) | divz G LfiVt, M)} endowed with the norm |u|r 2 := 

■^d 

r 3 ) + l‘ii^^lL 2 (Q ]r)) 3. Notice that M^) is a closed subspace of L|. 

Corollary 3.4. The space of admissible weak initial eonditions Aei is a closed subset o/Ljj^(ll, M^), 
and the space of admissible strong initial eonditions Ae 2 is a elosed subset of M^). 

Recalling the orthogonal projection 11, in (|^, we have: 
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Lemma 3.5. The norms \u \^2 and (\Ilu\‘jj + \u-n\‘^ _i j, are equivalent in U.^). 

^ (r, K) ^ 

Proof. Writing L^(r2, M^) = H Q) each v G L^(r2, M^) can be rewritten as 


V = Hi; + VP\;v, with / P\!V dfi = 0; 


( 22 ) 


in this way the mapping v i— )• P^/v G P[^{Q, M) is well dehned. Further, if divi; G M) we can 

see that P^/v solves the system AP^/v = divii, VP\7V ■ n = v ■ n. It follows that if i; G M^), 

then |VPv^^li 2 (Q^R 3 ) = (t” • n, Pyv)^_i^^ since P^/v is zero averaged in Q, there is a 

constant C > 0 such that IVPv^^li2(f^_R3) < C\v ■ ^^\'^Pvv\L^n,R3y, so |VPv?^lL2(n,E3) < 

□ 


C\v n|^_i < Ci|i;|j;^ 2 ,^(-q^ 


From now, for convenience, we suppose .4.=^ and endowed with the norm inherited from 

rd,v(S!. 1 ■ ■ " ‘ ~ 

Hilbert spaces. 


and from P[^^^{Q,, M^), respectively; from Corollary 


3.4 


the spaces A^i and .4=2 are 


Remark 3.6. In Section 2.3 we have considered the spaces of admissible conditions endowed with 
a suitable range norm also making them Hilbert spaces. Changing the norms now to those inher¬ 
ited from M^) and M^) will not cause any trouble concerning continuity properties. 

Indeed, we have that 4.=^ and Ae .2 endowed with the range norm are continuously embedded in 
M^) and M^), respectively, (cf. [Rodl51 Section 2.4]). Thus, from the completeness 

of both norms they are necessary equivalent (cf. [Brelll Corollary 2.8]). 


Next, we dehne the space 


Af:= 


z G 


d2M 


z = x^oPy± V(zi7rin + ZM+iTi) = 0 


M 

E 

2=1 


(23) 


Let Pj\f: —)• Af stand for the orthogonal projection in onto Af. Denoting the orthogonal 

subspace Af'^ to Af, we also denote Pj^± = /r 2 m — Pj\f: —)■ Af'^ the orthogonal projection in 

onto Af-^. Recall also the projections '■ i—;■ and i—;• onto the hrst M 

coordinates and onto the last M coordinates, respectively. We have the following property whose 
proof is given in the Appendix, Section |A.I[ 

PQ = QP for all (P, Q) G {Pm. Pm^} x {Qf. Qf}- (24) 


Notice that, for any u G .4=^ there exists at least one z G such that {u ■ n)njp = z = 
XIE^ ('^jP^Trinj; it follows, using (24), that the mapping 


u z'^'^ := Pmj-QY z = Qf Pm±z (25) 

is continuous from 4.=i onto and (u-n)njp = Indeed, to check that the mapping 

is well defined, we set another vector w G such that {u ■ n)njp = EQ^w, then necessarily 

EQ^[w — z) = 0 which means that Pm^Q^ {w — z) =0. On the other hand the continuity of 












12 


S. S. Rodrigues 


the mapping n i—>■ ii • n and the fact that both In • n| _i and | 2 :“'"L 2 m are norms in the hnite 

(r, K) 

dimensional space = {ix • n | it G (and, so necessarily equivalent) give us 


< C\u\ 




(26) 


Lemma 3.7. The norms \u\i^ 2 ^ (r 2 ,R 3 ) and (|nit||^ + are equivalent in .4=^; further they 

make A^i a Hilbert space. 


Proof. By Corollary 3.4, As.i endowed with the norm inherited from Ljjy(n, M^) is a Hilbert space. 

and from the fact that lit • n 


3.5 


The equivalence follows from Lemma 
equivalent norms in 


'f 


1 and I I xtr are 
□ 


Remark 3.8. Notice that the space M defined in (23) is not necessarily trivial, that is, it may 
contain nonzero vectors. See the Example in Section A.2 in the Appendix. 

_ 1 

Corollary 3.9. A element u G L^j,^,(n, M^) can he defined by (Hu, it • n) G R x Ha_A (T, M). 

_ 1 

Proof. Given {h^, h?) G R x M), and two elements and in M^) such that 

Hii^ = h^ = Hit^ and • n = /i^ = • n, then for u '.= we find it = (1 — n)it = VP^u 

and u • n = 0, which implies AP\ru = 0 and VP’^u • n = 0. Therefore it = VP’^u = 0, that is, 
= u^. □ 

Corollary 3.10. A element u G As.^ can be defined by (Hii, 2 :“'") ^ H x Q^Af-^. 

Proof. The mapping E: Af-^ —)■ L^(r, M^) is injective and Hz"'" = (it • n)n Ir- □ 

Now, let us be given four nonnegative constants 

0 < a < 6, 0 < e, 0 < 6, 

and two functions (p, (p & C^{[a, b], M) such that supp((/ 9 ) / 0, ip{t) > e for all t G supp((^), and 
(p{t) = 0 for t £ [a, a + (5] U [6 — (5, 6 ]. Further let 1 ? G ^^(r, M) be a function such that supp(i?) C O 
and i?(x) > e for all x G Tc; see ([Tl). 

Given a Hilbert space X, we define the orthogonal projection in Lp‘{{a, b), X)\ 

M 


Piif •= ^ f / fiT)^n{T) dr") a„ 

71 = 1 ^ 


where the cr„, n G Nq, are the eigenfunctions of the Dirichlet Laplacian := dtdt in (a, b): {(T„(t) := 
( 5^)5 sin(n 7 r(| 5 ^)) | n G Nq}, and {A^ = -(^)^ I n G Nq}; Atan = Kcin- 

Recalling the notations from Section 2.2 and inspired by an Example in |Rodl5l Section 5], we 
consider the auxiliary “control” space 

Qm ■■= p>x^^P^^PZPiiT^G\{a, b), r)|^ 

:= {Cl C = PX^oPx^PZPMiP'^vlo) and r? G G^((a, b), T)}. 


and consider the operator 


rj hG kZv ■= TX^oPZ±PmPm{t^v\o)- 


tO dO pO TDt 


(27) 
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Recalling the space if tv and the orthogonal projection IItv : H H^, see Section 2.2, from |Rodl5l 
Section 5] we know the following controllability result; 

Theorem 3.11. Let us be given u G and N € N, then there exists an integer M = 

1^1 ^ G No with the following property: for every vq G H, we can find rj = r]{vo) G 

G‘^{{a, b), r), depending linearly on vq, such that the control Q = K^rj = (px^'o P^±PMPM{‘f''^v\o) 
drives the system (16) to a vector v{b) G V such that IlNv{b) = 0. Moreover, there exists a constant 
Cuf, 1 , depending on |'u|vv;{a.b)|st, ip, ip, and b — a, but not on the pair (N, vq), such that 

l^(^o)lG2((a,b),r) ^ (28) 


3.2. Proof of Theorem 3.2, Let us fix a sufficiently large iV > 1 and let M be the integer given 
in Theorem 3.11 We organize the proof into main steps. 

(s) Step 1: driving the system, from v{0) = vq at time t = 0, to a vector i;(l) = vi ^ V at time 
t = 1. Let G Pjg-xQY'R^^ be the vector defined as in (25), and let cf) G ^^([0, 1], M) be a 
function taking the value 1 in a neighborhood [0, (5) of t = 0, and the value 0 in a neighborhood 
(1 — (5, 1] of i = 1, with S < ^. Let us also be given G . Then, the function K(/,(uo, = 

. . . , Hr,M+l^ M+2> • ' ' > «^r,2M) is in C^([0, 1], 


Next we consider the system ([16[) in (0, 1) x Ll, and the control f 


dtv + B{u)v — uAv + Vpv = 0, 


V 


p — 


div V = 0, 
i;(0) = uo; 


since ((up — EiEnp) ■ n) n|p = (uq • n)n|p — = 0, we have vq — EiPKp G H. By The- 

2 _L 1^- ,|2 


orem 
C 


2.5 


there exists a weak solution v satisfying the estimate 
[M,|n|^(o,i)|wE] + K\mao,i),R 2 M)), from which we can derive 








W I K.,-IJJ2A/ 


(29) 


Further, u|p vanishes in a neighborhood of t = 1, which implies that vi := i;(l) G H. Furthermore, 
since u|(q e yp(o,i)|st^ 

we actually have u(l) € V C Pf. 


from Lemma 


2.1 


(s) Step 2: driving the system from v{n) = G R at time t = n G Nq to a vector v{n + l) = Vn+i G V 
at time t = n+1, with \vn+i\‘jj < e~^\vn\‘jj- Now we consider the system (16) in (n, n+1) x Ll, and the 
control C = where = vi'^n) is given in Theorem 


3.11 


with (a, b) = (n, n + 1): 


dtv + B{u)v — vAv + Vp„ = 0, 




div u = 0, 

v{n) = Vn- 


Now, we observe that sin(m7r(t — n)) = (—1)”*"'sin(m7rt) and {\/2 sin(m7rt) | m G Nq} C ffQ((n, n + 
1), M) is an orthonormal basis in L?{{n, n + 1), M). For time t G (n, n + 1), the control 
can be rewritten as 

M M 

EE( rii,yn(^m{t)P^^'KiTa. + ri\„,a,n{t)Ti^ 

i=\ m=l 
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where rjf^ and are constants, and (Tm{t) = \/2sin(m7rt). Define : M. 2 MxM —^ n + 

1), M^)), mapping a matrix A = with real entries Aj^m £ M for j = 1, 2, , 2M and 

m = 1, 2, ..., M, to 

M M 

K^^^A := (/PxE? EE( 

i=l m=l 

and consider the space of matrices 


^fx ■■= 


A £ M2MxM 


= 0 


(30) 


we suppose M. 2 MxM ~ endowed with the scalar product {A, B)j^ := Ylm=ii-^i,rnBi^rn + 

^M+i,mBM+i,m)- It follows that A I—)• |P^xA|_yv( defines a norm in the range K^G‘^{{n, n + 

1), r) = K'^'^^M.2MxM- Since is linear and continuous from r/ £ G‘^{{a, b), F) into G‘^^{{a, b), F) 
(cf. [Rodl5l Proposition 5.1]), from the finite dimensionality of K^G‘^{{n, n + 1), F), it follows that 
\Pj^^A\j^ < C'|r 7 |Q 2 ((„ p), where A is any matrix satisfying K^^^A = K^rf. Moreover the map¬ 

ping rj A^ := Pj^±A from G‘^{{n, n -|- 1), F) into Af^ is well defined, that is, A^ is the unique 
element in J\f^ that solves := Kfr/. 

As a consequence of ([2^ we have that Defining, for each 1 < 

j < 2M, the functions K^{n, t) := X]m=i m•) = ^“(n, ■){vn) ■= 
{Ri{n, •), K 2 {n, •), ..., K 2 m(^i ')) is in G^{[n, u -|- 1], M^^); furthermore 






‘(n. 


and 




(n, •)('^'n)lcl([n,n+l], 


< <^riri 


] bnln- 


(31) 


Now, Lemma 2.8 and the continuity of the mapping v i—)> v{n -1- 1), from the space iy((n, n + 
1), M^), L^(D, M^)) into M^), imply that 


|u(n -h 1)1 






on the other hand from the definition of , in (27), we have that i;|p vanishes in a neighborhood of 
t = re+1 and, since K“(n, •) satisfies (31), we have that |u(n+l)|y < Cji^l ^ . i lunllr. Now we use 


, n + 1) |st 


yVnln- 


the fact that Ff 7 vu(n-|- 1 ) = 0 to o 
the Stokes operator (see Section 

then, for big enough N, such that > e^G 


itain a!Ar|u(n-|-l)|j^ < |i;(n-M)|^, where oat is the A^th eigenvalue of 




2.2), which allow us to write |i;(n-|-1)177 < a^/Cr,., i 

r, -, 1, we have that Vn+i ■= vin -|- 1) satisfies 

[l^lyy(n, n + l)|st J 

\'^n+l\\[ E e ^\Vn\\[- (32) 


@ Step 3: concatenation; a stabilizing control. First of all, we fix the functions (p and (p (appearing in 
Theorem 3.11) for the interval (a, b) = (1, 2) and then set (p{t) := ip{t — n + l) and (p{t) := ip{t — n + l) 
for t £ {n, n + 1). Since u £ (cf. (I!l))> the integer N in Step may be taken the same in 
each interval (n, n -|- 1), n G No; then, the same holds for the integer M in Theorem 3.11, with 
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(a, b) = (n, n + 1). Now we show that, given (vq, G x the control 

= r^(v K^\ ■= f if t G [0, 1); 

“ 0 ) tJ ■ EK^(n, ■)(vn), if t G [n, n + 1), with n G No; 


(33) 


stabilizes system (16) to the zero solution. Here, for n G Nq, Vn ■= v{n) where v is the solution of 
the system (16) in (0, n) x Q with control C*l(o n)xr- 

From ( |29| ), and the inequality 1 < for t G [0, 1], we obtain 

l^(^)lL2iv(n,R3) < + |k°|]^ 2 m^ , for all t G [0, 1]; (34) 

on the other hand for t > 1 we also have, |u(f )|^2 where [tJ > 1 denotes 

the biggest integer that is smaller than t, defined by 


[rj G Z and [rj + 1 > r > [rj, for all r G M. 


(35) 


Cji^l ^je''^e“^(^“^)|u(l)||^. Using ( p4| ) (with t = 1), we can conclude that 


Thus we obtain |u(t )|^2 jja^ < 


(I^o|£, 2 .^(q^r 3 ) + |«^rlR 2 A./j , for all t G [0, +CX)). (36) 

(s) Step 4: control estimate. Defining the mapping 


_ *,A/„ ^05 K?), if t G [0, 1); 

^ { R^{n, -Xvn), iffG[n, n- 


+ 1), with n G No; 


(37) 


(38) 


we see that the control in (33) can be rewritten as 

= Ek^’\vo, «°). 

Notice that for any giveiipositive integer n G Nq, the control function vanishes in a neighborhood 
of n. Indeed, from Step[^ vanishes in [1 — S, 1], and from Step[^ it also vanishes in [n, n + (5] U 
[n + 1 - 5, n + 1], because supp(<y9|j^^^_^^j) C supp((^|[^ ,j_,_;^]) C [n + <5, n + 1 - <5]. Let us be given 


-)MTn2M : 


into 


A G [0, A); then the mapping uo e-)■ e2*K“’''(uo, k)-) is linear and continuous, from A^i x Qj 
H^(Kq, M^^). Indeed, the linearity follows essentially from the linearity of system (16) and from the 


linearity of the mappings (vq, k,^) i— >• k^{vq, k°) and Vn i—t ^“(n, •){vn)- The boundedness follows by 
direct computations: we find 


e2 K 


ii,X 


e^'K^{vo, 




+ ^ e^K“(n, ■){vn) 


nSNo 


Hi((n,n+l),R2Af) 


< 


^[A] I '^T)|i^l((0,l),R2AT) + ^ 6 

nSNo 


Xn 


K^{n, ■){Vn) 


m{{n,n+l),R‘2M) 


^ ^[A,|ii|^st] ( l^olL2;^(n,R3) + |k°|k2AL + Y1 

neNo 


pAni |2 

e \Vn\H 
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and, using (36) (with t = n) and the identity XlneNo 


dA-A)n ^ _^A-a) 


62’k*,A (I'D, 


1 — 

2 


it follows that 


.0|2 


j^1(Ko,M 2M) - <^[(A-A)-hA,|,i|^.t] (l^oU2,^(n,R3) + K\R2MI, 
for any given 0 < A < A, which finishes the proof of Theorem 


3.2 


Corollary 3.12. The solution v = v{vo, k^) in Theorem 3.2 satisfies the estimate 

2 _ 

< C 


A 

e2 'n 




-i(r2,K3)) “ hA,|ii|y^,st] V^0lL|i^(r2, 


+ k 


0|2 


(39) 

□ 

(40) 


Proof. Proceeding as in the proof of Lemma 2.8, we start by noticing that w ■= e^'v solves sys¬ 
tem dlT] ), in each interval of time (a, b) C Mq) with {w{a), g, Kg) = ^e 2 “n(a) 

From Theorem 2.5, we have 


\ A A - . 

—^e2 ’n, r.e2'K’ 


\w\ 


W(in, n+1), Hi. (n, M3), H-l (O, 


52 , A. 


4- A|2 




< C„,, ., (e"”|„(n)|). 

thus, from (39) and (36) and from the continuity of the inclusion L^(n, M^) C 


w 




[a-a’'^’I“Iw=‘] 


< Cr 1 , 1 I ^e^"^e / e^*e ^*dt + l) (yf^,]R3) + |«^rlM2M ; , 


VnEN 




which implies (40), because X)neN + /mo dt +1 = +1, since A < A. □ 


4. The Oseen-Stokes system: feedback stabilizing control 
In this Section, we show that the finite-dimensional exponentially stabilizing control (cf. Theo¬ 


rem 


3.2) can be chosen in feedback form. In order to be more precise we will need to derive first 


some auxiliary results and consider a suitable extended “equivalent” system. 


4.1. Some auxiliary results. Once more we recall the projection 11 and the space Af, see 
and (23). 

Lemma 4.1. Let v solve system (16), with ( = Ek. Ife^'Ilv G L^(]Ro) H) and e^' k G Af'^), 

then e 2 'v G IF(Mo) M^)), with 


62 V\ 


iy(Mo,m. (o,m3),l/-i(o, 


-^[I«Iw=oA] (l^oli2.^(n,M3) + 1^2 -h |e2 


(41) 
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Proof. We rewrite 


A 

e 2 'v\ 


WiRi,H^.^{n, m 3), 77-1 (n, JJ3)) 


div 

from Theorem 


2.5 


as the sum le^ 
we can derive that 


+ 


a^ 2 i*^*i2 ii2 

^liv((0, l),77i (0,R3),77-i(n,M3)) <<^1^2 lci([0,1], R) TIm/((0, 1), 77]. (72, 




„ 77-1(72, R3)) 

2 


and, from Lemma |2.81 we have that for all t > 1 

which allow us to obtain 


di,(72,R3) ^ ^[liilwst] 1 )Il2.^(72,R3) + l«^l77i((t-l,2),Ar-L) 


(42) 


^'Il2(Ri,77];^(72, 
+ 00 

n=l 


+ 00 

n=l 
/•n+l 


L2((n,n+l),77i (72, 


^)lL2iv(0,R3) + I^l77l((2-1,2),A/'^) 


dt 


n=l 


+ 00 / .^+1 

2A / l„|(t-l) 


e24* Pv{t - (^^^3) dt + e^(’*+^)|7c|i 


77i((n-l,n+l),A^-L) 


- ^[I«lw=i]' 


-.2A I 


r+oo 


|e^*^Wli2^^(72,R3)dt 


+ 00 




n=l 


Since ez 'dfK = dt{e^' k) — ^e^'K, we can derive that 


+00 


^Il2(Ri, 771;^(72,R3)) < C'[|li|^,t,A] 


^[l“lwsl ’ 


e2 v\ 


L2(Ro,L2.^(72, 


+ ^ |e2 


((?^—1, n+l), A/*^) 


n=l 


— • 12 I — • I 

'*^lL2(Rg,7,2,^(n,R3)) + 2|e2 k|j^1(Rj,,_V^) 


Thus, from dt{e 2 'v) = ^e 2 'n + e^'dtv, and since v solves system (16), we can obtain the estimate 

which allow us to derive 


I — • 12 ~P^ 

I®" ^llV(Mi,77i;^(72,M3),77-1(72,R3)) - ^[|n|wst. 


e 2 v\ 


L^Ro,Ll.in, 


+ e2 "k 


77i(Mo,Ar-L) 


and then, using (42), we arrive to 


^+(Ro, 771 (72, R3), 77-1(72, 


^ '^[|ii|^st,A]i^l'^0ll2.^(O,]R3) + |e 2 7l|l,2(]Rg,,f^2.^(f^,]R3)) + |e 2 k\^hPRo,M^) 


4-,n2 
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Finally, from K{t) G , we have K{t) = Pj\f±K{t) = from Lemma 


3.7 


it follows that 

□ 


f(t )|^2 JI 3 ) < C{\P-'>j{t)\H + l<3/ which allow us to derive ( |4l| ). 

X Q, with ( = Ek and v{s) = Vg- If 


Corollary 4.2. Let s > 0 and let v solve system (16), in 
ez’IIi; G L‘^(S.s, H) and ea'ic G then 

|2 


A 

e2 ‘n| 








7Il2 (o,R 3) + I®^ + |e 2 ^I_H'1(R3,A/'-L) 


(43) 


Proof. Since v solve system (16), in x Q, we have that w = w{r) := v(r + s) solves system (16), 
in Mq X with w{0) = Vg, C = := EK{r + s), and Ug{r) := u{r + s) in the place of u. 

Since Irtslw®* < by Lemma 

|2 I |„G-t|2 


4.1 


we have that |e 2 , 

I iVK(Ro,//div(^d 

A 




is bounded by 


(l^-lL2^^(n,R3) + j, that is 


e2 


(-^)u|2 


iy(R^,Hi;^{o,R3),/f-i(o, 


^ <^[|iiUst,A] (l^slL2.^(n,R3) + 


L2(R„,H) 


+ e2 






which implies (43). 


(44) 


□ 


We will also need the following corollary of Theorem 3.2 and Corollary 3.12 


Corollary 4.3. Let us be given s > 0, A > 0 and u G W®*. Then there exists M = Cji^j ^ > 1 with 

the following property: for each {vg, k^) G Alsi x , there exists a “control” vector function 

i^u,x _ g L7^(Ms, m 2 M^ weak solution v of system (16) in x 11, with 

( = Ek^’^, satisfies the inequality 


A 

e2 ’u 


Vy(Rs,H^i„(n,R3),H-i(n,R3)) 
2 


+ |e2‘K“’^(u,,, 0| 


div 


Hi(Rs,R2Af) 


Vg 




_|_ |k.s |2 
I l^rlR2M 


(45) 


Moreover, the mapping (vq, k°) ^ k^’^{vo, k^) is linear. 

Proof. As in the proof of Corollary |4.2| we change the time variable t = r + s and can redu ce the 


problem to the cylinder Mq x H. W e may take (A, 2A) in the place of (A, A) in Theorem 
in the place of A in Corollary 


3.12 


3.2 


and A 


From Theorem 3.2 and Corollary 3.12 there is a linear mapping {vg, k^) i—)■ f,{vs, k^) = f,{vg, K^){r), 
r G Mo, such that the solution w = w{r) for system (16), with t(;|p = and with Ug{r) = u{s + r) 
in the place of u, satisfies {w, Q('^0(0) = ('^s, '^t) 


A 

e^'w 


U^(Ro,//i. (n,R3),H-i(Q, 


+ |et-e(u„ <)|^i(K„_ 


R2M) 


- (|^^Il2.^(0,R3) + K\l 
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Then we can conclude that {v, K){t) := {w, ^){t — s) solves system (16) in 

2 


(vs, K®) and that the estimate 


Cm 


) -^1 I'^slyyst] 


V |2 
Us r2 


e2 
s |2 


K—)u +|et(-*) 


IV (Ms, (n , R3 ) ^ H-1 , RT ) 


(o,]R 3 ) + l«^rlR 2 Mj holds, which implies (|4^. 


X VL with (v, Qi k){s) = 

l^{Vs-, «^r)|//l(]R^^]R2M) < 

□ 


4.2. The extended system. In order to be able to use the dynamical programming principle, we 
will need to rewrite system (16) in a suitable way. We start by observing that the mapping 

Fa-H u ^ ijlu, 

is continuous and surjective. From Corollary |3.10| it is also injective. Then by the Inverse Mapping 
Theorem (cf. [Brelll Section 2.3, Corollary 2.7]) it has a continuous inverse G £,{H x —?■ 

Msi). Further, since EQ^Af-^ C M)n we have that Q^Af-^) C M^), because 

any u G Q^Af^) C satisfies u = VP^/u, AP\ru = 0, and n • VP\/u = n • u = n • G 

775(11, M), which implies P\/u G H^{Q, M) (e.g., see |Tay97i Chapter 5, Proposition 7.7]). 

Now we rewrite system (16), in x ^ with = Ek, in the extended form 


dtv + B{u)v — kAv + S/py 

= 0, 

dtK = 

X, 

(46a) 

divv 

= 0, 

u|r = 

Ek, 

(46b) 

v{s) 

II 

, k(s) = 


(46c) 


where {v^, Kg, x) will be taken in 77 x A/"-*- x L^(Ms, A/"-*-). 


Remark 4.4. Notice that the compatibility condition (u(s) • n)n = EQ¥k{s), which is required to 
guarantee the existence of weak solutions for the “equivalent” system (16), is indeed guaranteed by 
the condition v{s) = F^^{v^, Q^Ks). Notice also that if we impose k{s) at the boundary, then we 
can impose only a tangential initial condition nu(s) G 77 for v{s). Essentially, the initial “weak” 
condition for the extended system (46) is the pair {v^, Kg). 


Observe that {v^, Kg) i—>■ Kg), Q^Kgj is an isomorphism in £(77 x Af^ —)• Msi x 

Qf^Al^). Form Corollary 


4.3 


it makes sense to consider the following problem: 


Problem 4.5. Let us be given s > 0, A > 0, u G W®*, and let M £ N be given by Corollary £3 , 
Then for given {v^, Kg) G 77 x Af'^, find the minimum of the functional 


P{Rs,F-^) 



X 

2 

X 

2 

A 

K, x) := 

6 2 TIu 

+ 

L 2 (Rs,H) 

e 2 'k 

+ 

L 2 (R,,Ar-L) 

02 


on the set of functions 

■= \^iv, K, x) £ Zy e 5 -( 


V, K, x) £ Z, 


s’^ and (v, K, x) solves (46a)-(46b)| 


satisfying A{v, k, x) := (nu(s), k{s)) = {v^ , Kg); where 

Zl := 1T(M„ m 3), M^)) x 773(M„ Af^) x L^{^g, N^) 

;= {(u, K, x) £ zl 1 v{s) £ MhJ. 
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Lemma 4.6. Problem 4-5 has a unique minimizer {v*, k*, x*). Moreover, there exists a continuous 


linear self-adjoint operator G C{H x Af-^) such that 

Ks))hxN^ = ^s) 


\R: 


A, s I 


l£(HxAri) < 


lie 


^As 


(47a) 

(47b) 


Proof. We suppose endowed with the norm inherited from and we start by observing that 
:= {(u, K, x) G I A{v, K, x) = {v^, Ks) and M^{v, k, x) < c} is bounded, for 
any {vf, Ks, c) £ P[ x Mq- Indeed, since the initial condition {v^ , Kg) and initial time t = s are 


the set A,,h ^ 


fixed, the boundedness follows from Corollary |4.2[ On the other hand, from Corollary 4.3 the set 
AvH,ks •— ^ I ^s)} is nonempty for any given {v^, Kg) £ H x Af'^, 

that is, the mapping A: Xg’^ —)• P[ xAf'^ is surjective. Further, we observe that M^{v, k, x) induces 
a scalar product in Zg’^: 

((u, K, x), {u, r], O) 


Mf 

( X A \ / A A 

:= e2 TIu, e2 ’IIu + e2'K, e2'r7 

V J P{Rs,V) ^ 


+ ( e 2 'x, e 2 '^^ 


L2(R„,Ari) 


We can derive that Problem 4.5 has a unique minimizer {v*, k*, x*) = (v*, k*, x*)(v^, Kg), which 
linearly depends on (v^, Kg) (cf. |Rodl5l Appendix, Lemma A.14 and Remark A.15]). 

Again from Corollary |4.3[ we have that the mapping 


{iVs'\ (wf’^ ^s)) ^ («, ^s), (v*s, K*, X*){v^’‘^, At^))^A 

is a symmetric continuous bilinear form on 77 x Af'^ which is bounded by ^ ^ on the 


As 

unit ball; thus, the optimal cost can be written as (47a), where R£ is a bounded and self-adjoint 
operator, which norm satisfy (47b). □ 


Next we consider another minimization problem related to Problem 4.5 


Problem 4.7. Let us he given s > sq > 0, A > 0, u G W®*, and let M £fi he given by Corollary U-^ 
Given (u^, Kg^) £ H x Af'^, find the minimum of the functional 

2 


N. 


S0,S\ 



X 

2 

X 

2 

A 

K, x) := 

6 2 'IIu 

+ 

LR{so,s),H)) 

02 'aC 

+ 

L2((S0,S)W-L) 

02 'X 




+ ((nu(s), k{s)) , (nu(s) k{s))^ 


on the set of functions 

{v, K, x) £ X := |(u, K, x) £ 2^g^ g-j I V solves (46a)-(46b)| 
that satisfy A{v, k, x) := (nu(so), k(so)) = {v^o^ Kg^); where ■= 


Reasoning as in the proof of Lemma |4.6[ we can derive that Problem 4.7 has a unique minimizer 
(^*o,s) ^so,s){'^?o’ ^«o)) which is a linear function of (u^, Ksq) £ H xAf~^. The following Lemma 

is the dynamic programming principle for Problem 4.5 (with s = sq). 
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Lemma 4.8. The minimizers of Problems \4-5\ and \4-T\ have the following properties: the restriction 
o/«o’ to the time interval (sq, s) does coincide with 


^50/? 


and the restriction of {v*^, x*sq){Vsq, Ksq) to the half-line = {s, +cx)) does coincide with 
(vt, Kt, x:)(ni;:^(s), Ki^is)). 

Proof. From „ Kso) ^ ^ and {v*, k*, ^(s), k*^^,{s)) € , we neces¬ 

sarily have {v*^^g{t)-n)n = EQ^K*^ g{t) for any f G [sq, s], and {v* (r) - 11)11 = EQ^K*{r) for any r > s. 
Thus, since ^(s) = we find that ^(s) • n = v*{s) ■ n, and from ^(s) = nu*(s), it 

follows that = Xj\Uv*^^gis), Therefore, the 

concatenation 

/ «,s, <o,s: '^so){t), if t G [so, s]; 

is a function in . Analogously, we can see that also the concatenation 

'^*so^ ^*so)i'ofo^ if t G [so, s]; 

'^*soi^))it), if t G> s; 


(T, K, >c){t) := 


is a function in xl^. 


By the definition of (u* , k* , x*) and (u*, k*, x*) we can conclude that 


Msoiv*,, <, >c*so){vZ, Kso) < K, x) = <, ^^o)l(.o,s) 




so 


so’ 

2 


Lmso,s),H)) 


+ 




e2 K 


so 


H 


I,2((so,s),Ar-L) 


+ 


-t . 

ePx, 


so’ “so’ 
2 


so 


L2((so,s),Ar-L) 


+ ' (n<(s), <(s)) , (n<(s) <(s))) < (48) 


4.5 


it follows {v*, K* Ka)(vZ, Kso) = (v, K, x). 


From the uniqueness of the minimizer for Problem .............. v'^so’ 

In particular, k*^, «^so)I]r, = «, ^s)(n<(s), k*o(s)). 

On the other side, from the definition of (n'o^s, «^*o,s, ^*o,s)('^'i^, ^so) and (48), we also have 

= ^so,siv*so,s^ <o,s’ Ko,s)iTso^ /^so) < <o’ ^so)(^^’ '^So)l(,o,s) 

= ^4^0^^’ = ^s^o(^so’ ^so’ ^so)(^i^’ «^so) < K, x). 

Necessarily, it follows that (T, k, X) = x*^){v^^, and Kso)I(so,s) = 

(^S0,S, ^S0,S’ ■^S0,s)(^S 0’ ^So)- ^ 

4.3. Linear feedback stabilization for the extended system. In this section we prove the 
following Theorem 4.9 which says that if we see x as our control and (u, k) as the state in (46), then 
X can be taken in linear feedback form. 

Theorem 4.9. Given u G and X > 0, let M = Cy\u\y^,\] G N 6e the integer constructed in 
Theorem 


3.2. Then there is a family of operators G C{H x J\f-^ —>• such that the following 

properties hold: 

(i). The function s 1 —)• s G [0, + 00 ), is continuous in the weak operator topology, and it 


holds 


K 


A, s 
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(ii). For any given sq > 0 and Ksq) ^ H x the solution of the system (46) with x = 

X Q, and satisfies the estimate 


/C 


A, 11 


(nii(t), K{t)) exists, in 
2 

(v, k) 
for all a> sq. 




(7 


|(nu(a), K(a)| 


HxjV-^ 


(49) 


Proof. We organize the proof into main steps. In Step we use a Lagrange multiplier approach to 

of Problem 


derive two key optimality conditions for the minimizer ^so,s) of Problem 4.7 In Step|2| 

we use those conditions and the dynamic programming principle to hnd the linear feedback rule. In 
Step 1^ we show the uniqueness of the solution under the feedback controller, prove the bound of the 
feedback operator norm, for each instant of time, and prove estimate (49). Finally in Stepj^we prove 
the continuity of the time-dependent family of feedback operators in the weak operator topology. 


(s) Step 1: Karush-Kuhn-Tucker Theorem. First we notice that, considering as usual H as a pivot 
space, we can extend the projection n: —)• H to a mapping II: —)• W by simply 

setting (n/, u)v',v •= (/) '^)all tt G F (cf. beginning of Proof of Theorem 5.3 
in |Rodl5| ). Then, we define the spaces 

Aio := VF((so, s), ^ ^ 

Ts := {(u, k) I (v, K, 0) G T’o, u(so) G Ae^, and v{so) = T'_^^(ni;(s), Qfnis))} , 

X-.= Xs X L\{so, 

y-.= HxM^x L\{so, s), V') X L\{so, s), M^) x GI^q{{so, s), F), 


where we denote 


G'iv,o((so, s), L) := {u G Gav((so, s), F) | u(so) • n = 0} 


(where we may understand u(so) ■ n as (Eiu}(so) ■ n, with the extension Ei as in Section 2.1). 


Next we define the affine operator F: X ^ y,hy 

F{v, K, x) := (nu(so) - k(so) - Ksoj ^{dtv + B{u)v - nAv), dtK - x, u|p - Ek). 

We show now that the derivative dF of F, dF{v, k, x) = F{v, k, x) + (u^, 0, 0, 0), is surjective. 

Indeed let us be given (yi, y2, 2/3, 2/4, 2/5) £ V- Setting K{t) = 2/2) for all t G (sq, s), x = —2/4, and 
taking K in system to be the inclusion from g((so, s), F) into Gav(('S0) s), F): g C = Krj = 
n. then bv |Rodl 5 l Theorem 2.11], we have that svstem (©, in (sq, s) X n, has a weak solution v for 
the data (uq, g, C) = Qfy2), 2/3, 2/5 + ^2/2), because necessarily {Efi^{yi, Qfy2) - Efiy^ + 

Hy2)(so)) • n = -{ EQfy2 ) • n = 0. 


Therefore, we find dF(u, k, x) = ( 2 / 1 , 2 / 2 , 2/3, 2/4, 2 / 5 ) and can conclude that dF is surjective 

in Problem 


4.7 


is a minimizer in X, that is 
the Karush-Kuhn-Tucker The- 


We see that the minimizer of the cost 
iVso,s, J ^ and satishes F{vl^ „ „ x*^ J = 0. 

orem (e.g., see [BRSlll Theorem A.l]), there exists a Lagrange multiplier 

(/is, 0s, Qs, TOs, 7s) G y', 

with y = H X X L^((so, s), V) x L^((so, s), x q((so, s), F)', such that 


<o,^’ ^so,s) + (Ts, Os, Qs, ros, Js)odE{v*^ „ k* 


so,S’ "^SO, s. 


= 0 . 
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Hence, for all (z, <;) G and all ^ G L^((so, s), Af-^) we have 

0 = 2 f e^\Uvl^^„Uz)Hdt + 2 f dt 

JSQ JSQ 


(50) 


HxN^ 


+ ((n^(so), ?(so)), (/^S, Gs))hxN^ + / df 


/so 


+ / (n(5t2: + - zzAz), gs)v/y df + (7s, 2;|r - ^1./ 

./sn i^o,r^-^o,i 


0 = 2 


/so 




(>^so,si OaT-l dt+ {Ws, -OaT-l dt , 


(51) 


'so 


where for simplicity we denote Gq p := o((®o, s), F). 

(s) Step 2: properties of the optimal triple. Letting z run over all z G VL((so, s), V, V), with 
2;(so) = z{s) = 0, and taking ? = 0, then from (50) and from M^) = H'©{Vp | p G L^(H, M)} 

(see, e.g., [TemOll Chapter 1, Section 1.4]), we can see that for some pg^ G L^((so, s), L^(H, M)), 

- dtqs - vA.qa + B*{u)qs + Vpg^ + 2e^^Ilv*^ ,,{t) = 0 (52) 

where B*{u) is the formal adjoint to B{u)'. dehned for given q G V and v G M^) by 

{v, B*{u)q)L2(^n,R3) := {B{u)v, q)H-^n,R3),H^,{n,R3 y Furth er, we suppose that we have hxed an 
appropriate choice for the pressure function pg^ (cf. |Rodl5l Section 3.2]). To fix ideas, let us choose 
Pg^ satisfying f^pg^ dH = 0. 

From G L'^{{so, s), H) and qg G L‘^{{sq, s),V), it follows dtqs £ .^^((soi •s), HOj ™ particular 

qs G G([so, s], H) (cf. |LM72( Chapter 1, Theorem 3.1]). Using again (50), with arbitrary z G 
hF((so, s), V, V) C G([so, s], H) and ? = 0, we derive that gs('So) = hs and 

qs{s) = -2R^;J(nu;,,(s), <„,,(s)) G H 

where, recalling that i?^’^(nu*Q ^(■s), ^so,s(s)) G H x Af-^, we define 


(53) 


<?(n<,s(s), ^:o,sis)), Rt%{uv:^,sis), K:^,sis))} ■■= Rt’^{nv:^,s{s), k:^,s{s)); 


->A, s 


A, s / 


(54) 


'Rl’liUvl^^sis), ^:,,s(^)), <1(nu:„,s(s), '^:o,s(^))j €HxAf^. 

Since Husjj s £ L‘^{{so, s), H), by standard arguments we can prove that gs £ hF((so! s —2e), D(L), H) 
for any 0 < e < taking, again in (50), arbitrary 2 supported in [sq + e, s — e] with 2 ; G 

^((^ 0 ) s), M^), M^)), and ? = 0, we obtain, using (52), 

0 = ((n • V)qs-Pg,n, ^lr)L2((so,s),L2(r,R3)) + (7s, z|r)(GB = )/,Gj;- (55) 


On the other hand, relation (51) implies that 

Ws = 2e^*x* 

which shows us that zus is inc 
[so, s]), because by Lemma 


4.8 


(56) 

ependent of s (or, more precisely Ws{t) is independent of s, for t G 

H 


we have k^o) = (^^ 0 ’ < 0 ^ '^so)l(,o,,) 
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Now taking z = 0 and ? G Hq{{so, s), AA-*-), from (50) we obtain 

rs rs 

0 = 2/ dt + / (9t?)_v-Ldt +(7^, 

Jso Jso 

and, introducing the adjoint E* : (Gg’p)' —)• J\f-^ of S: J\f-^ —)• G'g’p, defined by 

(H* 7 , k)j^± := ( 7 , Ek)f^i,ss, ^i,s , for all k G AA"*-, 
'.'-"o.r-’ ’'-' 0 , r 


we arrive to 


“*7^ = - dtws = 2e^*(l - - 2Ae''*x*g_^. 


Xt/ 


.At 


(57) 


Then, we see that % is independent of s. Since the space {/ G Gg p I supp(/) C (sg, s)} D {/ G 
G°°([sg, s], -f7^(r, M^)) I supp(/) C (so, s)} is dense in L^((sg, s), T^(r, M^)), by using ( [5^ we 
obtain that (n • V}qs — Pg^n = js is independent of s. Here, the last equality is to be understood 
in 

Suppose now that we are given another si with sg < s < si. Then, the difference (q, pg) := 
(qs - , Pqs - Pqs-, ) satisfies 

-dtq- uAq + B*{u)q + Vpq = 0 

in (sg, s) X H. From the observability inequality in |Rodl5l Inequality (3.4)], we obtain 

\q{s - 2e)\‘jj < C\pqn - (n • V)g|^i((,_ 2 ,_,)^p), = 0 for all 0 < e < 

Since q G G([sg, s], H), it follows that g = 0 in [sg, sj. In particular, we can conclude that qs 
and Vpq^ do not depend on s and that 


- 2R^’*(nu*j,_^(s), K*s^^s(s)) = qs{s) = qsXs) £ for all s G [sg, si] 


(58) 


which, in turn, implies that qs G VF((sg, s), D(L), H) and pq^ G L^((sg, s), F7^(H, M)) (cf. |Rodl5l 
Section 3.2]). In particular, if we choose pq^ satisfying J^Pqs dT = 0 then the identity 7 ^ = {n-V)qs — 

Pq^n is meaningful in the subset L?‘((sq, s), Hi^(T, M^)) C (Ggp)'. 


Now, we observe that the function 77 := i^*sq,s + ((' “ '®)^*o,s(^o) “ (• “ 'So)^*o,s('S)) is in 

Hq{(so, s), AA-*-) and, from (57), solves the Poisson equation 

(1 - At)Tq = \e-^-E*-is + Ax*^ ^ + ^ ((• - s)k*o,s('®o) - (• - so)k*s^, sis)) , (59) 

with At := dtdt- Since the right hand side is in L^((sg, s), J\f-^), by standard arguments (for el¬ 
liptic equations) it follows that rj G 77Q((sg, s), AA-*-) n 77^((sg, s), AA-*-), which implies that ^ G 
H'^iiso, s), Af-^). 

Whence, taking z = 0 and an arbitrary ? G H^{(so, s), AA-*-) in (50), with ?(sg) = 0, and using (54) 

- 


and 


we can arrive to 0 = 2 R^’ 2 (nu* ^(s), k* sis)) + zusis). Thus, from (56), we obtain 


>^so,sis) = -e ^"R^ 2 (nu*^^^(s), K*o,s(s)), and Lemma 


4.8 


>^Lis) = (s), K* (s)), with JC^’" := -e 


leads to 


A, s ._ —\g dA, s 


u, 2 ' 


(60) 
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(s) Step 3: uniqueness, operator norm, and exponential decay estimates. We know that given 
Ksq) ^ H X Af-^ the pair (v, k) = (v*^, k*^)(v^, k^q) solves (46), in x with the feed¬ 
back control ( |60| ): 

dtK = /C^’' (Hr;, k), (61a) 

f Ip = Ek, (61b) 

K(so)=Ksg, (61c) 

This solution is unique: indeed, if (v, k) is another solution with (Hu, k)(so) = (v^, Kgg), then 
(d"^, d'^) = (v — V, K — k) solves 

dtd^ + B(u)d^ - nAd^ + Vpd- = 0, divd’' = 0, d^'lp = Sd^, 

(d^(so), d-(so)) = ( 0 , 0 ), d- = f IC^^’^iUd'^ir), d-ir))dr, 

JSQ 


dtv + B{u)v — vAv + Vpv = 0, 
divu = 0, 

viso) = Qfl^so), 


and, from (47b), ( [54| ), and ( [60| ), we have 

-A, r 


\ICriIlFir), d^r)r^, < I]|(nd’^(r), d^{r))\i,^^,. 


Then, for the function z{t) := JC^’^{Ild^{t), d'^{t)) 




, using Theorem 


2.5 


and d^{so) = 0, 


^ ^[Nw-.A,i] + \d^ir 


a] 


JC^’\Ud%s), d^is))ds 


'so 


= c 


[l'“lw=t ' a] 


'so 

pr 


+ (r - So) / z{t)dt 

H'-((0,r),jVA dso / 

2 \ 


/CT^(nd"(s), d"(s))ds 


'so 


+ (1 + r - so)2;(t) dt 


AT-L 


(62) 


-^[NwsoA,!] [(i-so) J^^zis)ds + {l + r-so)z{t)j dt. 

Hence, for any r > sq and r G [sq, r] we have 

and from the GronwalHBellman inequality (see |Pac73( Theorem 1]) we obtain z{r) = 0. Since r > so 


can be taken arbitrary, we can conclude that z = 0 which implies d'^ = 0. Using again Theorem 2.5 
we arrive to d’' = 0. 


The uniqueness of the solution of system (61), implies that we need to prove the estimate (49), 


for the optimal trajectory (u* , k* ){vg, Ksq) solving Problem 4.5 (with so in the role of s). In 


this case, from (47) and for any a > so, we already know that ) < 
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|(n<(a), «:o(«))lkAr^- Hence from (|44]), 


|et(-“)z;* |2 




_(_ I P 2 (' K* I ^ 

|e ^Solj^l(Ra,A^-L) 


- ^ ‘^*^^SolL2(Ea,W) U 1®^' '^’^dolffl(Md,ArJ-) 


!(-.)„. |2_ 


'^|K.iw..,Al e"<W)lL^A +e-"“(l + a2)m^(.,;i,. , , <uj) 

C|w,„A,iil(n";(“).'A;(“))ikAfx, 


that is, estimate (|4^ holds. 


@ Step 4: continuity in the weak operator topology. By definition, {/C^’* | s G [0, +oo)} C C{H x 
AA-*- —)• AA-*-) does depend continuously on s in the weak operator topology if for all {w, G {H x 
AA-*-) X AA-*- and sq G [0, +oo), ^)jg-± goes to Oa/"-*- ^ Soes to sq; with s G [0, +oo). 

Recalling (j^ and (60), we see that this property holds if it holds for the family {R^’^ \ s G 


[0, +oo)} C C{H X AA“*“ H X AA“*“). 

Given a pair {w^, w'^) G (i^ x A^-*-) x {H x AA-*-), we can write 2{R^’^w^, w‘^)jjy..j^± = {R^’^{w^ + 
w"^), ~ ~ we see that it suffices to prove 


that 


{Rl'''w,w)Hy^_^± as s So, for any w e H x Af-'-, 


?A,So 


(63) 


(64) 


or equivalently (cf. |Con85[ Chapter IX, Proposition 1.3.(e)]) that 

{Rl'"°w,w)Hycjvr- as n ^ + 00 , 

for any sequence (<5,i)„gN of real numbers, with 0 < < 1, <5n —>• 0, and any w £ H x AA-*- (still, with 

s = So ± (5n > 0). 

We consider separately two cases: s \ so and s Z' sq. To shorten a little the notation we denote 
zl^{w) := (IIu*^, kI^){w), for a given w £ H x Af-^. 

(a) The case s \ so- If s = sq + (5n, we write 


R£ w,w 


HxN^ 


(65) 


= (R 


;)A, SoH“< 5 t: 


{w - zl^{w){so + 5n)),w'^ ^ + (^R^’'*°+'^"z*g(u;)(so + <5n), wj 


HxN^ 


Rewriting the last term in equation ( [6^ as (^R^’^°~^^^z*g^{w){so + 5 n), w — z*^(w)(so + 6 n)^ 

''°^'^"Zso(^)('SO + <5n), 2l^iw){so + 6n)'j 


' HxN^ 


which, recalling Lemma 


gramming principle), can be written as 

+ <5n), w - zl^{w){so + hn)j 


4.i 


+ 

R'xA/'-L 
(the dynamical pro- 


HxM^ 


+ (w 


HxjV^ 
2 


|e 2 1®^ '^So(''^)lL2((5p^So+5n)W-L) 1^2 >^Soi'^)\L2((so,So+5n),J£-^) 
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Using the self-adjointness of ^e arrive to 




,iv] — (r^’^°w, vci) (66) 

= - zl^{w){so + Sn)),W + zl^{w){so + 

- \^^"^<o(^)\haso,so+S„),H) - <i^)\lH{so,so+S„),M^y 

Now, for the first term on the right-hand side of (66), we have 

^o+'^"(u, - zl^{w){so + 6n)),w + z*^{w){so + 


^[l^lwst, A, A]® 


HxA/'-L 

Hso+Sn)\^ _ z* {w){so + Sn)\HxjV^\'^ + + ^n)\HxM^ 


and, from the continuity of zl^{w){t) in the time variable t, we can conclude that the this term goes 
to zero with 5n- On the other hand, the remaining (last three) terms on the right hand side of (66) 
also go to zero with 6n, because the triple k*^,, >^sp)(ri^) do not depend on 6n- Therefore we can 
derive 

(67) 


lim {rI’^°^^-w,w 

n^+oo ' ^ 


HxjV-^ 




HxjV^ 


(b) The case s sq. Though we will follow the same idea, this case carries a few additional 
difficulties. If s = sq — (Jn, we write 

and, writing z* c^(t(;)(so) = z* c^(t(;)(so) — rc -|- tc in the last term, we obtain 


HxN^ 


SO — Sn 

R^ w,w 


HxA/'-L 




= \(^^'^V*so-SnH\hiiso-S^,so},H) (68) 


HxJ\f^ ' 

2 


+ |e2’Kso-<5„(w^)lL2((so-<5„,^o),Ari) + \^^'^*so-Sn('^)\h{iso-S„,so),M^) 
+ {Rl’"°iKo-sA'w)iso) - w), zl^_^Jw){so) - ivj 

+ 2 Zso-<5„(«^)(so) - V?j^ 


HxjV^ 


HxN^ 


Now byl^we have, in particular. 


I ^ 'N 12 I ^ \ 12 

I C([so —Sq —(5n + l]) I C'( [SQ — , SQ — (5n + l], ) 

^ '^[|ii|^st,A, A,so]l^lHxAr-L 
which implies, since 5,1 < 1, 

|e" ’n<0-5„ (^)li2((so-5„, so), H) + |e2 ■ 

^ <^[|ii|^st,A, A,so]l^l^xAri'^>2- 


(69) 


(70) 
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On the other hand, from x*^{w) = {Uv* .^{w), k* .(w)) and 


|e" ■ ^so-S^ MlL2((so-Sn, so),Ar^) 

^ (l^^"^^so-SnMlh((so-Sr,,so),H) + |e"'<o-< 5 n(^)lL 2 ((^o- 5 n,so),Ari) 


and so, from (68), (70), and © , we arrive to 


lim ^''‘w,w 

n^+oo 


HxjV^ 


(^^’^°(^so-S„HCso) + U’), z;^_^^(w)(so) - 


HxN^ 


( 71 ) 


(72) 


= w)^ 


HxM^ 

Denoting {w^, w^) := w £ H x A7-*-, since V is dense in H, there exists a sequence Wn = 
(tcjf, 'u/^)n£n, with £ V, such that \wn — w\hxJ£->- to zero as n goes to +oo. Thus, writing 
w = w — Wn + Wn and using (69) we obtain that 

Nso-5„(^)(so) - w^IhxAT-l 

< Ko-sA'^ - Wn){sQ) -w + WnlnxM^ + |Zso-5„(w^n)(so) “ WnlnxM^ 

- ^[l“lw=t.A, ~ '^nlnxjv^ + |Zs(,_5„(w;n)('So) “ WnluxjV^J 


and 


< 


^°(Zso-5„(«^)('So) + w), zl^_^Jw){so) - w'j 
^[|ii|^st,A, A,so]l'“^lHxAr-L (|a^ - U’nlffxAT-L + l^so-'Sn('^AAo) “ Wnlj/xArA 


HxN^ 


which implies that 


lim 

n^+oo 


bA^°{z*^_s^{w){so) + w), zI^_^^{w){sq) - u;) 


R'xA/'-L 


(73) 


< '^[liiU.t.A, ^lirn^ |z:o-5„(«'n)(so) - Wn\HxM^- 

It turns out that = Vn{-) '■= v*^_g^{wn)i-) — w^ solves ( [l7| ), in (sq — 6n, sq) x D, with g = 


B{u)wX — i^Awn, K = E, g = K*^_g^{wn){t), and Vn{so — 6n) = 0. From Theorem 
b(i)lH-i(n,K3) < C'(I^WIi-^(n,M3) + AAnlni (Q,K3), we obtain that 


2.5 


and 


Ko-5AWn){so) - wX\l2 


< C 


Vi2 


= l^n(so)li2 


[l“lw(»o-^n..o)|wk] \LH(so-S„,so),V) + Ko-sA^n)\Hiaso-S„,so),M^)) ; 

so 

Sq 


on the other hand we also have (?/;„)(•) — v/^ = ^ >t:*^_S^(wn)(t)dt and 

2 9 

|«^^0-<5„(^n)(so) - <^n \xl^_sAWn){t)\j^^ dt. 

J So <577, 
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lim 

n—I'+oo 






(74) 


Then, by ( fTO] ), ( [7l| ), and (73), we obtain 

+ w), z*^_ 5 ^(u;)(so) - wj 
^ A, i,so] <^n {\wX\v + 

= ^[|rU.,A, iso]l^lHxArx 

Recall that the inclusion R C 77 is continuous. Notice that will go to +oo with n, if G H\V. 

We need to show that the sequence {wX)n&n can be chosen such that 6n\wX\y —>■ 0, as n —>■ +oo. 

We know that the sequence Hn+iw^ G Hn+i C V converges to , in H, where IItv: H —)> 
is the orthogonal projection in H onto the space 77jv spanned by the first N eigenfunctions of the 
Stokes operator L (see ([Tol)) but, we have no guarantee that 5n\^n+iw\y —)■ 0, that is, Hn+iw^ may 
be not a good choice for w^. Next we show that this issue can be overcome by somehow “slowing 
down” the convergence of to , in H: we define the sequence iV: N —?• No, n i—?• by 


Nn ■■= 1, 

Nn ■= j, 


if 

if 


202 ^ 
j > 2 and 


< 


{j+l)aj+i ' — jaj 

where(aj)jgis}p is the nondecreasing sequence of (repeated) eigenvalues of the Stokes operator (cf. Sec¬ 
tion 2.2). We can see that Nn goes to -l-oo with n, because 6n goes to 0. Hence, setting Wn '■= HTv^rc^, 
we have that —)• , in H, and \wX\y < C\'VwX\‘j^ 2 (^Q r9) = C{LwX, wX)h < l^n Ilf 

then 6n\wX,\y < C6nOiNn\wX.\’H < for all Nn > 2; thus 5n\w^\y —)■ 0 as n —>■ +oo. 

From (72) and (74), we arrive to 


lim 

n^+oo 


t^A, So — 

w^w 




= ( R^’''°W, w 


HxjV-^ 


(75) 


Finally, from (67) and (75), we see that (64) holds. That is, the family depends continuously 


on s, in the weak operator topology. 
This ends the proof of Theorem |4.9[ 


□ 


4.4. Miscellaneous remarks. As we have said in the Introduction, looking for feedback finite¬ 
dimensional controllers supported in a small subset of the boundary, is motivated by the importance 
of such controllers in applications. We give a few remarks concerning this point. 


4.4.1. Dimension of the controller. The range Em of the controller depends only on the norm lulyyat 
of the targeted solution u, and the feedback rule depends on time. We do not address here the 
problem of finding an estimate for M, that is of crucial importance for application purposes (e.g., 
numerical simulations). For internal controls, this problem has been started in |KR,15b( IKR,15a) 
in the simpler case of the ID Burgers system in a bounded interval (0, L), and estimates on the 
number M of needed controls is given that depend exponentially in M^ef := + u~^X )2 in 

the general case, and that are proportional to M^ef in the case of no constraint on the support of 
the control, with W = L°°(Ko, L), M)). However also in |KR15b| the results of numerical 

simulations suggest that it might be sufficient to take M proportional to Mref also in the general 
case. Following |KR15b) . we see that in the particular case, the estimate for M is derived from an 
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inequality like aM > where aM is the Mth eigenvalue of the Laplacian operator. Hence we 

may (perhaps) conjecture that in the dD case, d G {2, 3}, it should be enough to take a number M 
of internal controls proportional to because um ~ uCiMd (cf. |Ily09| ). Does 

the conjecture hold true? Can we derive similar estimates for the case of the boundary controls we 
treat here? These questions will be addressed in future works. 

Recall that, in the case of stationary ft, for example in |BT04[ IBTlll IRTIO) . we can find rather 
sharp estimates, though M does depend on it and not only in the norm lulwst- The method cannot 
be (at least not straightforwardly) used in the nonstationary case. 


4.4.2. Lyapunov functions. Once a feedback control is constructed, it is easy to find a time-dependent 
Lyapunov function for the problem in question. Indeed, the functional 

poo 

<P{r,w)= \TA'Sr,tw\1j^j^^dt 

J r 

decays along the trajectories of (61), where Sr^w denotes the solution {v*, Kl){t) of (61), in x 
with the initial condition (Hu*, K*){r) = w: we may write 

poo poo 

(HnS, kI){t)) = / dr = / | JU5o..(nuS, k5)(0)|^,^x dr 

Jr Jr 

poo 


from which, together with (49), we can obtain 
d 


ds 


<^(s,(nnS, K*)(s)) = -\{Uv*„ Ac5)(r)|^^^, 


< - 


s=r 
T=7 -1 


Y J- I ^ ^ ^ ^ ^ I ^ 


< -( 


Another Lyapunov function is the “cost to go” from time t = r onwards, that is \L{r, w) = 
{K^w, w)jjy^j^±, where R” := is the operator defining the optimal cost. Indeed, from the 

dynamical principle we have 


{Uv*„Kl)ir))=^iO, iUv*o,Kl)m 

X 2 

- ea'Huo 

^ L^{{0,r),H) 


e2 Ko 


L2((0,r),Atr) 


e^'dtKl 




which implies, using (47b) and 
d 


ds 


s=r 
Xr 


4>(«. (nt,;, «;)(»)) = 

< -e"|(n»j. 4){r), (no;, K;)(r))„„^o 

It is, however, difficult to write down the functions L' and L/ m. a. more explicit form. 
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4.4.3. Riccati equation. In the case of internal controls it was shown in [BRSlll Remark 3.11(b)] 
that the operator defining the optimal cost, and from which we can obtain the feedback law, satisfies 
a suitable differential Riccati equation. For applications (e.g., simulations) it is important to have 
such equation at our disposal; for example, we refer to [KR15b| where the optimal feedback rule has 
been obtained by solving (numerically) a similar differential Riccati equation, in the simpler case of 
the ID Burgers system and internal controls. 

It turns out that also in our case the operator R := i?* := from which we can obtain the 

boundary feedback rule, satisfies a differential Riccati equation. 


First of all we recall (54) and rewrite R = {Ri, R 2 ), where 

H), 




(R(u;^ w'^), {z^, z'^))HxAr^ =■■ [ 2 ^ z^] 


1- 

1- 

1_ 

R2,1 R2,2 



Ri e £(H ^ H), R2 €C{HxM^ 

and for any two pairs (tc^, w"^) and (z^, z^) in H x , we have 

{R{w^, w^), (z\ 2^))iyxAr-L = (Rl(w^^ w'^), z^)h + {R2{w^, w^), z‘^)j^±. 

Now observe that we can further decompose (Ri, R 2 ) and define Riq G £{H), Ri ^2 £ R(AA-*- — H), 
R 2 ,i £ £{H —AA-*-), and R 2,2 £ R(AA-*-) by 

{z^, Ri,iw^)h ■= ((^\ 0), R{w^, 0 ))hxM^ > Ri,2w‘^)h ■= {{z^, 0), R(0, w‘^))hxMz , 

(z^ R2,iu^^)At-L := ((0, z^), R{w^, ^))hxM^ > (^^ R2,2w‘^)mi- ■= ((0, z^), R(0, w‘^))hxMz ■ 

It follows that we can write 

n _ -^ 1,1 -^ 1,2 
“ [R2,1 R2,2. 

with z^Ri^jW^ := (z^, Rijw^)h and z‘^R 2 jw^ '■= (z^, R 2 ^jW^)j\fz. 

The operator R := e“^*R satisfies, for all t G Mq 

R - RA - A*R - RBB*R + C = 0 

with 

A := 

where stands for the Oseen-Stokes operator 

L^(i): m 3) ^ n(-izAu + B{u{t))v), 

and is defined by 

^ Ru{t)RQY 

Observe that, from Lemmas |4.6| and |4.8[ we have 

{R%Uv*„ kI){s), {Uv*„ k 5 )( s))^,^x = kI, dtKl) 

from which we can derive that, formally, 

{ds\s=tR"(^V^, K*Q){t), (nU(), K*Q)it))Hy,j^i_ 

+ (-R*5sL=t k*q){s), (nUf), K*Q){t))H^^i_ + {R\llv^, K*o){t), (Bu]), k*o){s))h^^± 


(76) 


Ru{t) ^u(t) 

A 

'i o' 

; B := 

'o o' 

; C := 

'1 o' 

0 0 

~ 2 

0 1 

1 

0 1 

") ^ 

0 1 


e2‘nuo(t) - e2*Ko(t) 


H 






AT-L 
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where = dtKl- From ds\s=t (Hfo, = {-Lu{t)vl{t), K{t)), recalling @ and using Lu{t)vl = 

QYKq) = Q^Kq), in the matrix form notation we can write 


[nus 

+ [Uv*, 






+ [lluo A^o 


-(LuT 


KqI R 


L{i 

.-A- 


-e ^'R 2 ,i 


^u(t) 

-e ^■-^ 2 , 2 . 



1- 

* c 

R 


L J 


-e-^'R: 


= -e^’ [nu^ / 

Therefore, we can conclude that 

0 0 

0 0 

-^2,1 R2,2 


1 

0 

1—1 


1— 

*0 

R 

0 1 


kI 


— e 


2,1 

* 

2,2 

-A 


R 

UvY 


L ^0 J 


[nus 


Kf, 


1 

* 

0 


0 

0 

_ 1 


*0 

R 

0 R 2 , 2 . 


-^ 2,1 R 2,2 


Kg _ 


R-R 


= e-^'R 



0 

* 

s 

A. 

1 o' 



Ft -\- G 

0 1 


+ e 


-A- 


d R 
0 R. 


2,1 

'* 

' 2 , 2 _ 


R — e 


-A- 


d i 
0 R 22 


0 

-^2,1 


0 

R2,2 


Using the symmetry of R, we have that R 21 = Ri ,2 and R 22 — ^ 2 , 2 , and we can derive 


e-^-R 


0 0 

^ 2,1 R2,2 


+ e 


-A- 


d R 
0 R. 


' 2,1 

I* 

' 2 , 2 _ 


R — e 


-A- 


d i?: 
0 R. 


2,1 
'2,2j 


0 0 

;R2,1 R2,2 


-A- 

0 Ri^2 

0 

0 

e 

0 R2,2 

^2,1 

R2,2 


= e-^'RBBR. 


Hence, from (77) and B = B*, it follows that 


R-R 




1 

* 

0 

0 0 


L(T(o)‘ »J 


R + e^'C- e-^'RBB*R = 0 


which implies that for R = e R, using R = e R — AR, we have the identity 


R-R 


Rib 

0 


u(t) 


{LuY 


R + C - e-^^'RBB*R + AR = 0, 


Hi,i 

R2,1 


R * 

2,1 

R2,2 


, the feedback control, in Theorem 


4.9 


(77) 


(cf. system (61), and (60) 


which is equivalent to (|76). 

Once we have R = 

) is then given by 

X = }CY ' (Hu, «)(•) = —R2,inu — R2,2«^- 

The structure of the feedback control is comparable with that proposed in [BadOOt below Equa¬ 
tion (1.22)], in the case of a stationary targeted solution. We also remark that here, taking advantage 
of the finite-dimensionality of the controls, we consider a simple dynamics dtn = x on the boundary, 
while in [BadOO] Equation (1.12)] a different dynamics is proposed to deal with a larger class of initial 
conditions, but with controls that are not necessary finite-dimensional. 

Solving (76) numerically, say with finite element method as in [KB,15b] . can become a quite 
demanding problem as the number np of mesh points increase. In the case of internal controls we 
will end up in solving a matrix Up x Up-dimensional problem. On the other hand, in the boundary 
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case and considering the extended system instead, the problem will not get much worse because the 
dimension of the problem just increase by the number M of controls, and M will (or, is expected to 
be) be much smaller than n^. We may expect that the numerical demand in solving the resulting 
matrix (rip + M) x (rep + M)-dimensional problem is comparable with that in solving a matrix Up x Up- 
dimensional problem. 

Notice also that the “less standard” operator k i—)> Q^k), from J\f-^ to A^i can be con¬ 
structed by solving M = dim < M Laplace equations of the form 

= 0, n • Vps = n • Hit, 

(for each R G and then assigning k) = VpqM^. 

Concerning, the numerical solution of the Riccati equation we refer the reader, in particular, to 
the works in |BBSW15[ IBenOGl IBLM97] . See also |BK84[ iBSZOSl IKM90j . 


5. Stabilization of the Navier-Stokes system 


In this Section we prove that the feedback controller constructed in Section to stabilize the 
linear extended Oseen-Stokes system (46) to zero, also stabilizes locally the corresponding extension 
of the nonlinear system (15) to zero. The main result of the paper is given in Section [5.2[ before, due 
to some well known issues related with the existence and uniqueness of solutions we need to recall 
some definitions (cf. [Rod 141 Section 5]). 

5.1. Solutions for the nonlinear systems. Let a, 6 G M be two real numbers with 0 < a < 6. 

Definition 5.1. Given / G L^((a, 6), M^)), k G L^{{a, b), L®(n, M^)), and yo £ H, we say 

that y G b), R) n b), H), with dty G L^((a, b), V), is a weak solution for system 


dty + {y-^)y + B{k)y-i^Ay + Vpy +f = 0, divy = 0, 

y\r = 2 /(«) = yo, 

in (a, b) x H, if it is a weak solution in the classical sense of |Tem01l Chapter 3, Section 3.1]. 
For simplicity, and for r > 1, A: > 1, we define the subspace 

•= m 3)) n L”((a, 6), L'=(L!, R^)). 


(78) 


Definition 5.2. Given h G b), H ^(fl^ M^)), z G 0^)^ and uo G L^jy(n, M^), we say that u 

satisfying 

u G L^{{a, b), m3)) n L^{{a, b), m3)), dtu G L\{a, b), m3)), 

is a weak solution for system JT), in (a, b) x 11, with 7 -|- C = -^Ir '^(®) = if y = u — 2 is a 
weak solution for the system (78) with f = h + dtz -|- {z ■ V)z — vAz, k = z, and yo = uq — z{a) G H. 

Definition 5.3. Given u G b), M3)), 2 ; G acid vq G M3), we say that v 

satisfying 


u G L\{a, b), m3)) n L“((a, 6 ), dtv G L\{a, b), H-\n, m3)) 

is a weak solution for system (El), in (a, b) x H, with ( = zjp and v{a) = uq, if y = u — z is a weak 
solution for the system (78) with / = dtz + {z ■V)z — iyAz + B{u)z, k = u + z and yo = vq — z{a) G H. 
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Analogously, we define the global solutions in Mq x let us denote, for simplicity, Ai, 

{/ I f\{o,T) G A’((0, T)) for all T > 0}, where ^((0, T)) is a suitable space of functions defined in 
(0, T) C Mo- 

Definition 5.4. Given h G L^q^(Mo, M^)), G 0^’^ and uq G M^), we say that u 

is a weak solution for system in Mo X n, if n|^Q is a weak solution in (0, T) x fl, for all 

T > 0, for the same system with the data h | j-q , /c | 2 ^) • 

Definition 5.5. Given u G ((0, T), M^)), z G and vq G M^), we say 


that u is a weak solution for system (15), in Mo x fl, if u|^q is a weak solution in (0, T) x fl, for 
all T > 0, for the same system with the data ^|(o t) ’ ^l(o T)- 

The existence of the solution in Definitions 


|5.2| and 5.3 can be proven following classical 
arguments, and some continuity observations as mentioned in [Rodl41 Remarks 5.1 and 5.2]. Then 


it also follow the existence of the solutions in Definitions 5.4 and |5.5[ Concerning the uniqueness, 
following an argument as in the proof of Lemma 5.3 in |Rodl4j . we have the following: 


Lemma 5.6. The solution in Definitions 
and 


5.1 


in Definitions 5-4 


5.2 and 

5.5, is unique if it is in 


5.3, is unique if it is in The solution 


, loc * 


5.2. Main Theorem. To have enough regularity to deal with the nonlinear problems, we take 


strong admissible initial conditions in Ae2 = {u £ 


u\ 


= Ez for some z G A/’"*'} C Aei 


Recall that Ae 2 is the set of admissible strong initial conditions for the Oseen-Stokes system (see 
Section 2.3). 

Now we rewrite system (15), in Mq x D with f = Hk, in the extended form 


dtv + B{u)v + {v ■ V)v — uAv + Vpv = 0, 

divu = 0, 

u(0) = Qfno), 


dtK = fC^’fiUv, k), 


V r = 


k ( 0 ) = Ko, 


(79a) 

(79b) 

(79c) 


where {v^ , Kg) will be taken in HxAf-^ x L^(Ms, and is as in Theorem 4.9 (cf. system (HI))- 

Theorem 5.7. For given u G W®* and A > 0, let M = ^ he as in Theorem 

exists e > 0 with the following property: if 


4 . 9 . Then there 


{vq, ko) G (iLniL]i^(D, M^)) X A/"-^, with F/{vq , Qf ko)\j. = Ekq, 


(80) 


and 


Qfno) 


M, 




< e, then there exists a weak solution {v, k) in the product space 


VLioc(Mo, iLj;.^,(D, M^), L^(D, M^)) x AA-*-) for system (79), which is unique and satisfies the 

inequality 

t>Q. (81) 








Before the proof we need some auxiliary results. 

Lemma 5.8. If {vq , kq) satisfies (80), the solution v in Theorem 4-9 satisfies 


A -,.|2 




r>0 
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Proof. We know that, taking a = 0 in (49), 
2 


e^' {v, k) 


< C 


K-o) 


HxjV^ ’ 


(82) 


since w := e^'v solves ( [l7| with g = —^v, g = e ^'k = e^' + Jq ^^’^(nf(s), K(s))ds^, and 

K = 'E, from Theorem 12.61 we have that 


w 


W({r, r+l),Hj.^{n, R3), L'^{n, 


^[l“lws‘] 


w{r)\ 




+ \v\ 




On the other hand, by constru ctio n, in the proof of Theorem 4.9 we know that the solution coincides 

with s = 0; and, from ([^ it follows that \v\l2^(^r,r+i),L^. (n.RS)) ^ 


with the minimizer of Problem 

|2 


4.5 


A,, 4] ^o)\hxM^- Therefore, using 9*7? = + e 2 -9*^, 


\^\w{{r,r+l),Hl.jn,R3),L2{n, 

- + \ivo , ^o)\Hx^f^ +Mq{v, k, dtn)^ ; 


and from Lemma 14.61 we have 


\w 


W{{r,r+l),Hl^{n 


,R3),L2(0,R3)) - ‘^[|il|w=t,A, 4] (l^(0l^i.^(n,R3) + I 




In particular, |7n(LrJ + 1 )|^i^^(^^k 3) ^ ij 

[rj G N as in (35), and 


PPl 




+ |(^^^> f^o)\H^^±), with 


I 12 I 12 

l^llV((r,r+l),/72.^(n,R3),L2(f7_]R3)) < \w\w{ilr\, [rJ+2), (H, R3), L2(f7^ 

- ^[l«lw=tO, x] (l^(L^-l)lHii^(n,R3) + l^(M + l)l^i;^(n,R3) + 2 > ^o)\Hx^f^ 

- ^[|71|w=t,A,4] (e^^'’V(H)l^i.^(o,R3) + '4o)|^xAri) • 


If [rJ >0, that is, if r > 1 by Lemma 2.8 we also have 


^HIh1;^(0,R3) < {\v{[r\ l)li2.^(Q,]R3) + l«l/7l(([rJ-l,LrJ)WV) 


(83) 


(84) 


and from 7? = e ^'k and 9*7? = ^t? + e 2 ■9*7t, we can conclude that 




+ l^i^li2((LrJ-l, PDW-L) 
< + |e-t*(9*7? - i??)li 2 ((L,j_i, Lr-JiW^) 


( 85 ) 















36 


S. S. Rodrigues 


thus, from (|83|), (|84|), and (85), and recalling (47a) and (47b), it follows that if r > 1, then 
|2 


A 

e 2 u 


iy((r,r+l),H2 (Q,R3),L2(Q, 


= \w\ 


Ty((r,r+l),//2 (Q,R3),L2(Q, 


< 1] (e^W|^(L^j _ + IK^, «o)|LAr^ 

Now, ^ implies |u([rj - i]e-^(W-i) ko)|L^x, thus 

W((r,r+l),//2.^(0,R3),L2(Q_]R3)) < , A, i] ’ '^o)|HxAr-L ’ for r > 1. 


e 2 V 


From (83) we also have 


•„.|2 


e2 V 


M/((r,r+l),//2,^(0,R3),L2(Q_R3)) < , A, (l^olHii^(0,R3) + > ^OT r e [0, 1). 


Since {vn , kq) satisfies (l80|), we have IkoIat^ < Ci |Hko| i < C 2 |i'(0)|rri /q ]r 3 '| and 

— div'^ ’ > 

|e 2 ^^l^((r,^+i),H2.^(o,]R3)^i2(f^^K3)) ^ I'^( 0 )|^i.^(f^_jj 3 ) , for T > 0 . 

which ends the proof. 


□ 


Inspired by Theorem 4.9 and Lemma 5.8, we define the Banach space 

:= 1^ G (Moj ^^div(^) kl.z^ < oo| 

endowed with the norm |z|^a := sup^>o 


A 

e2 'z 


lV((r,r+l),R'2.^(0,R3),L2(0, 


For given constant /9 > 0 and {vq , kq) satisfying (80), we define the subset 

■= ^ I ^( 0 ) = Qf Ko) and \z\^x < p| 2 ( 0 )|^i.^(f^_K 3 )} > 

and the mapping T: Z^ —>■ Z^^^, z z, where (z, k) solves 

dtz + B{u)z - izAz + Vpz^z = -n((7 • V)z), dtK = K,^' {liz, k), 

div z = 0, ^\r ~ 

2;(0) = , Qfno), k( 0) = Ko, 

Lemma 5.9. Under the hypotheses of Theorem B there exists p > 0 such that the following 
property holds: for any 7 G ( 0 , 1 ) one can find a constant e = > 0 such that, for any {vq , kq) 

satisfying (80) and \iF^^{vQ , (o.rs) < the mapping T takes the set Z^ into itself and 

satisfies the inequality 


( 86 a) 

( 86 b) 

( 86 c) 


\'^{zi)-'^{Z 2 )\z^ <l\zi-Z 2 \z>- for all zi, Z 2 G Z^ 


(87) 


Proof. We divide the proof into main steps: 

@ Step 1: a preliminary estimate. Consider the system 
dtz + B{u)z — izAz + Vpz = f, 
div 2 = 0 , 

z(o) = Qfno), 


dtK = )C^’'{Uz, k), ( 88 a) 


z|p = r.K, 
fi:(0) = Kq, 


( 88 b) 

( 88 c) 
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where / G H). If {z, k) is the solution of system (88) with / = 0, by Lemma 


5.1 




(89) 


We are going to derive a version of this estimate for suitable nonzero /; for that we denote by 
Sq ^{vq , kq) the solution (z, k) of (88). In the case / = 0, the operator Sq ^ is linear; by the 
Duhamel formula we can write 


{z, K){t) = Ko) = Ko) + / S'^tifis), 0)ds 

Jo 


(90) 


where S{, Kg) denotes the solution of the system (88), with the initial time moved to t = s, 
and the initial condition {v^, Kg)- On the other hand, from (49), it follows in particular that 

then 


'«)(i)li2.^(n,M3)xA7r < /«o)|V2 / 5°j(/(s), 0)ds 

^ 0 




[llilyyst, A, 


\{vo, ^^o)\HxM^ + |(/(s), 0)|^^^T ds^ 


(91) 


Now we can find, again with [tj G N as in (35), 


j-t W ^fc+i 

/ e2^|/(s)|/i-ds < V / e“2^e^®|/(s)|Hds 

Jo Jk 


hj 

sE 

k=0 


< sup 


rk+l 


As 


ds 


rk+l 


e^^1/(s)lHds 


jeN 


A w 

A^"\f{s)\Hds] ^ 


k=0 


e-^" ds 


and for the sum of the series, we can find X)fc=( 


k+l^-Xs 


< Er=o 




-^k _ (l-e-Qs _ 


^+1As , 
k 

t 


ds = 


C J1j SUPjgN 


o<i<hJ 


A2(l-e“^) 
1 
2 


= Ojij. Hence we obtain the inequality e2^|/(s)|j:^ ds < 


A^^\f{s)\H dsj and, recalling (91), 


/ 


^Xt 


\{z, K)(i)|^3^xA7-L < C'[|tl|^,t,A, 


lAo, /«o)| 


2 

J7xA7-L 


ri+1 


V 


+ sup 
i6N 

o<i<hj 


.,2As 


l/(s)l^,,ds (92) 
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for all i > 0. Now we use Lemma |2.8| to obtain 
|2 




(93) 


+ (7 


[l^lw®*] 


K(r) + / /C^’^(nz(s), fi;(s))ds 




Recalling the bound in (62), we find 


k(’’+ l)lHii,(n,M3) < [\4r)\As, + \f\LH{r,r+l),H) 


+ ^[Nw-.A,i] l(n^i ^)WlHxAr^- 

LI IW - ’ AJ 


From (92), and (94), we obtain 


(94) 




/ 


^[FIw=* ’x] 


lie 


Ar 


nfc+l 


\ivo , i^o)\hx^^ + sup / e "|/(s)|^ds|; 

fceN Jk 

0<k<yr+l\ 


because the inequalities |/li2((^_^+i)^ < e e^®|/(s)|^ ds < e e^^^\f{s)\'jj ds < 

2e~^^ supfcgpj ^‘^^^\f{s)\‘jj ds do hold true. 

0<fc<Lr+lJ 


For t G (0, 1), from Theorem |2.6[ (62), and (92) we can also obtain 


(96) 


+ C 


[Flwst] 


k( 0) + / /C^''(n2;(s), k(s)) ds 


10 


//1((0, l),Ar-L) 

- ^[Flw=t] (|^(°)Ih1j^(RK 3) + l'^(0)liri + I/Il2((0,i),h)^ 


which, together with (95) give for all t > 0: 


/ 


\zit)\ 




lie 


—Xt 


rk+1 


in-R3)+ sup / e2^"|/(s)|^ds 

div'- ’ '' L.z-'FM I 7 . 


\ 


ken Jk 
0<fc<pj 


where we have used 0 < e = e^e '^(’’+1)^ and 1 < e^e for r > 0 and t G (0, 1) and also the 
fact that |K(0)|_yy^± < C'|u(0)|j|^i j^s) because {vq , kq) satishes (80). 
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Using again Theorem 2.6 and proceeding as above, and using \ K{t)\j^± < C 
now derive 


l^lTV((r,r+l),/f2.^(Q,R3),z,2(Q_]R3) < (II R3) + \f\L2{{r,r+l),H) 


+ C 




K{r) + / /C^’^(n 2 ;(s), k(s)) ds 


m{{r, r+l),Ar-L) 


we can 


x] 


lie 


—\r 


rk+1 


^(0)lHij^(n,R3) + ^^sup 

o<fc< W 


,2As 


\fis)\Hds ; 


which implies, since dt{e^'z) = ^e 2 'z + e 2 'dtz, that 

|2 


sup |e 2 5;(i)|pi/((r,r+l),H2.^(n,R3),L2(r2,I 


^ ( I^(°)I/^L(0,R3) +sup 


(97) 


rk+l 


keN Jk 


'\f(s)\Hds 


which is a version of (89) for nonzero /. 

(s) Step 2: 'I' maps Zp into itself, if |2;(0)|j|^i is small. We will replace / by —n((z • V)z) 

in (97). First we recall some standard estimates for the nonlinear term; from the Agmon inequality 

'- _ 1 1 

(see |Tem97t Chapter II, Section 1.4]), |u|ioo(f^^K3) < C'ikl|fi(o,R 3 )kl|^ 2 (Q^]R 3 ) < C'2|u|^2(j^_r 3), we 
can obtain \ {z ■ V)h)|i2(f2,M3) < C\z\h 2 (o,r 3)Iwllm. (o,R3) and 


div'' 
A , 


sup / e^^'Kz-V)z|^ 2 m K3)ds < sup |e^"2(s)|^i ^ ^3) sup 

fceN Jk se[o, pj+il ’ k&, 

0<fc<pJ 0<fc<pJ 


rk+l 




ds. 


Thus, inequality (97) with / = —n(( 2 ; • V)z) gives us 




li z ^ Zp, then 


^;;M2 



(98) 

Hi;^(Q,R3))N(0)l^ii^(n,R3) 

(99) 


and if we set p = 2C and e < -, where C = Cj|^| ij is the constant in (99), then we obtain 
C(1 + p^e^) < p if |2(0)|j^i, ]r 3) < e, which means that T(z) G Zp. 


@ Step 3: 'h is a contraction, if\z{0)\fji (f2^R3) is smaller. It remains to prove (87). Let us take two 
functions ^ 1,^2 £ Zp and let (T(zi), ki) and ('I'(z 2 ), '^ 2 ) be the corresponding solutions for ( [8^ . Set 
e = ^1 —^2 and {d^, d'^) = ('I'(zi) —'I'(z 2 ), ki —K 2 ). Then {d^, d'^) solves (88) with (d^, (i'')(0) = (0, 0) 
and / = n((z 2 • V)z 2 ) — n((zi • V)zi). Therefore, by inequality (97), we have 

_ rt+l 

|4/(^-l) - 'L(^- 2 )||a < ^^[|.|^.„A,i] sup / 


e^^^\{z 2 ■ V)Z 2 - {zi ■ V)2;i||2(Q^]R3)ds. 


( 100 ) 
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Straightforward computations give us 

\{Z 2 ■ V)z 2 - {zi • V)zi||2(f^_]R3) = |-(e • V)z 2 - {zi • V)e|^2(Q^]R3) 

and, from ( |100| ), it follows |^'(zi) - ^'(z 2 )||a < A, i] (l^il|^ + I^2 ||a) |e||A, that is, 

\'^{zi) — '^{z2)\2,\ — ^j2p|z(0)|^l^^^j^^]g3^|zi — Z 2 I 2 A (101) 

Choosing e > 0 (smaller than the one chosen in Step and) such that 2C2pe^ < 7^, where C 2 = 
^[\u\ t A A] constant in (101), we see that if |-z(0)|j|^i (f2^]K3) < e, then (87) holds. 

The proof of Lemma 5.9 is complete. □ 


Proof of Theorem \5. '1\ Form Lemma 5.9 and the contraction mapping principle (cf. |Smi83t Chap¬ 
ter 10, Theorem 6.1]) it follows that if u(0) = ^ I'^o) £ ^^) satysfies (80) and is 

sufficiently small, |'y(0)|j:^i (q_r 3) < e, then there exists a unique fixed point v G Zp for T. It follows 
from the definitions of T and Zp that together with a suitable function k, the pair (u, k) solves the 
system (86), with z = v. 

Since we can write {v ■ V)v = n((u • V)v) + Vpa, for a suitable function G 77^(11, M), it follows 
that if we set py := py^y + Pv then we can conclude that (v, n) solves (79). 

From ([^, taking z = v = z, we find |t>||A < a, A](l + (o,r 3)’ which im¬ 

plies (81). 

It remains to prove the uniqueness of (u, k) in the space 

f := VFioc(Mo, R^)) x 77iL(Ko, D Z\ 


Let {v, k) € Z he another solution for ( |79| ), and set (z, 9) := {v — v, n — k). Then z solves (15) with 

■u -|- u in the place of n, z{h) = 0 and ( = E9 = E /C^’^(IIz(s), Q^9(s)) ds = E Fj\^z{s) ds. 

3 ^ 

We can extend a given function g G 77av(F, M^) to the solution Fg ^ M^) of the Stokes 

system 

-/\Fg + Vpg = {),d:wFg = {),Fg\^=g (102) 

and the mapping g 1—)• Fg is continuous (cf. [TemOll Chapter 1, Proposition 2.3]). Then, we can 
define the extension F: EH^^^fRo, Af-^) —)• H^^^fRo, M^)) by 

M 

FEK{t) := Ki{t)FxRo P^rTTiU + KM+i{t)Fx^oTi. 

i=l 


Notice that F is injective because Fg 


is finite-dimensional, we have 


FEK{t) 


Moreover we can see that y := z 
0, and / = dtFE9 + {FEO • V)FEe 


= 0 implies g = Fg\-p = 0. Observe also that, since F: EAf-^ 




< C 


FEK{t) 


As 


— FE6 solves system ([TS] ) with k = u + v + FE6, y{0) = 

— 12 AFE 9 + B{u + v)FE6. Hence, by standard arguments 
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and ^ 


(y, F~e) 


HxAs 


= 2{dty, y)H + 2{FEdtO, FE9 )_a^^ we can find 


df 

<Ci 


(y, FEB) 


HxAs, 


U + V + FE0 


L°°(n,R3) 


IvI'h + 


-/ + dtFE9 


v 


+ 2{FEdt9, -y + FE9U^^ 


and from 


/ - dtFEB 

2 

< C 2 

FEB 

4 

+ C 2 

FEB 


V' 


As, 





+ C2 Iff + v\j^a 


'(n,R3) 


FE9 


{FEdt9, -y + FE9)_A^^ = {FEF^FFA{y + FE9), -y + FE9)_A^^ < C 3 (y, FEB) 


A^, 

2 

HxA-e 


we arrive to 

dt 


(y, FEB) 


HxA'= 


< Ci { Ifi + ul 




+ 


Since {v, B) ^ Z we have that g = C 4 + v\j^oo(^q j^s) + 

from the Gronwall lemma it follows that 

2 




FEB 


+ 1 


(y, FEB) 


HxAs, 


As, 


+ 1 ) is locally integrable and 


(y, FEB)it) 


HxAs, 


< e 


Jo ff(*) ds 


(y, FE9){0) 


HxAs, 


= 0 . 


We can conclude that v = v + z = v + y + FEB = v, and also that k = k + 9 = k, because 
FE:M^ ^ FEM^ is injective. That is, the solution {v, k) for (|79[) is unique in Z. □ 


5.3. Integral feedback rule and stabilization to trajectories. Given vq G A'E 2 C 
there exists a unique kq G Af-^ such that uolp = Let (u, k) be the solution of (79), with 

{vq , QyKo) = (IIuo, QYkq) = F_avo. Then we observe that v solves ( [I^ with ( = Ek and u(0) = 
T'_^^(nuo, QYKq) = Vq. From Theorem 5.7, we can conclude that if |uo|j:^i is small enough, 

then the integral feedback rule 


i;|p = Hko + H/ JC^’^Uv, QYK){s)ds = vo\^+ E I /C^^JUu(s)ds 

Jo Jo 

does exponentially stabilize system (15) to zero with rate That is, setting kV := fYJ^A 
the feedback control Q = {uq — h(0))|p + E JJ K^’®(tt(s) — -u(s))ds does exponentially stabilize sys¬ 
tem ([^ to the targeted trajectory u with the same rate provided uq — u(0) G As 2 F M^) 

and |uo — ^( 0 )|j:/i. (fj ]r 3 ) is small enough. 

Furthermore the corresponding solution u for system ([^ is unique in the affine subspace u + 
Wioc(Mo, m 3 ), m 3)). 


-A.s - 
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A. Appendix 


A.l. Proof of Property (24). From the definition of S in (14), we have that for any k G 

Ek = 0 if, and only if, [EQ^k, EQf^ k) = ( 0 , 0 ), 

because EQ^k G A^(r, M)n and EQf ^k G -b^(r, TF) are orthogonal in Lp'iT, K 
K G AA if, and only if, {Q^ k, k) G AA x Af, which implies that 

QPx = PmQPm, for Q G {Qf, Qf^}. 

Hence, since for any given (^, k) G x we have 

{PpQPjV^^l k)]r2M = {Pj^±^, QPmI^)k.‘ 2M = {Pj\f±^, Pj\fQPj\fK)^2M = 0, 
it follows that PmQPjV-'- = 0 ^-rid 


p2M 


This means that 


PpQ = PpQPp + PnQPm^ = QPp-: 

Pm^ Q = Q - PpQ = Q - QPp = QPj\f ^) 


for all QG{gf, Qi^}, 


which is the statement of (24). 


□ 


A. 2. An Example concerning Remark 3.8 We illustrate the fact that the space N defined 
in (23) is not necessarily trivial. ^ ^ 

Let O = (0, vr) x T^; the family {vTj | i G Nq} contains the family {cj,i | 
n G No}, with fTn(r, s) := :^sin(nr), (r, s) G (0, vr) x ~ (0, vr) x [0, 27r). 

Define the indicator operator TT 27r J : C{0, M) —)• LP‘{0, M), sending / to 

' TT 27r 'N 




f 22 l^/ defined by ^\f(r, s) := 

{’.32 V .3 ’ .3 2 


f{r,s) ifrG(f, f) 


0 


ifr G (0, 7r)\[|, -f] 



Now, setting the mapping x •= P(^ 2 l (Sus — ug), from direct computations we obtain the i 


iden¬ 


tities = 0, = X(|^ ( 3(73 - ag), and = X^^^ ^ (dua - cjg); we can check that 

X G C‘^{0, M) and supp(x) = [|, x C O. Now we can show that the functions xIE^H^(cj 3 n) 
and xIE^P^(crgn) are linearly dependent: we find 


d^X — 


92 


A’fx(o' 3 n) = 0-311- 


( 0 - 311 , XJi)l 2 (c>,r) 


Xn = 0-311 


10 


xn, 


(fTgn, XtP)l‘^(o,R) , 1 

-xn = o-gn -h Yq An, 


Ixnp 


2 Xx(o' 9 n) = o-gn- 

"'■■'LpO,R) 

from which it follows 3P^ ( 0 - 311 ) — P^(o- 9 n) = 30-311 — ^xn — o-gn — j^xn = 0. Therefore, if the 

functions (o-an) and (o-gn) are in the family {xIE^P*^(Tijn) | i < M}, then the 

XX X 

family is linearly dependent; it follows that Q^M C J\f contains nonzero vectors. 
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